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Abstract 

The purpose of this article is to introduce and study a relativistic motion whose 
acceleration, in proper time, is given by a white noise. We deal with general relativity, 
and consider more closely the problem of the asymptotic behaviour of paths in the 
Schwarzschild geometry example. 
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Corpora cum deorsum rectum per inane feruntur, ponderibus propriis incerto tempore ferme 
incertisque locis spatio depellere paulum, tantum quod momen mutatum dicere possis. Quod nisi 
declinare solerent, omnia deorsum, imbris uti guttae, caderent per inane profundum, nec foret 
offensus natus, nec plaga creata principiis : ita nil umquam natura creasset. 

Lucrecius 

1 Introduction 

The classical theory of Brownian motion is not compatible with relativity, as it appears 
clearly from the fact that the heat flow propagates instantaneously to infinity. A Lorentz 
invariant generalized Laplacian was defined by Dudley (cf [Dl]) on the tangent bundle of 
the Minkowski space, and it was shown that there is no other adequate definition than this 
one 1 , as long as Lorentz invariance is assumed. An intuitive description of the associated 
diffusion (i.e. continuous Markov process) is that boosts are continuously applied in random 
directions of space. We show that this process is induced by a left invariant Brownian 
motion on the Poincare group. The asymptotic behaviour of the paths of this process was 
studied (cf [D3]). 

Considering the importance of heat kernels in Riemannian geometry and the extensive 
use that is made of their probabilistic repesentation via sample paths, it is somewhat 
surprising that Dudley's first studies were not pursued and extended to the general context, 
namely to Lorentz manifolds. It is indeed easy to check that the "relativistic diffusion" can 
be defined on any Lorentz manifold using a development, as done below. The infinitesimal 
generator is the generator of the geodesic flow perturbed by the vertical Laplacian. But 
such an extension would have little appeal, if some natural questions such as the asymptotic 
behaviour and the nature of harmonic functions could not be solved in some exemples of 
interest . 

Here we provide a rather complete study of this question in the case of Schwarzschild 
and Kruskal-Szekeres manifolds, which are used in physics to represent "black holes". 
The specific interest of these manifolds comes from the vanishing of Ricci curvature, their 
symmetry, and the integrability of the geodesic flow. 

The picture that comes out in the Kruskal-Szekeres case appears quite remarkable, with 
paths confined in a neighborhood of the singularity, while their velocity increases, and an 
infinity of SO3 -invariant harmonic functions. 

One difficulty of the study (and it might explain why Dudley had few followers) is that 
no explicit solution was found. The reason is that, even after reduction using the symme- 
tries, the operator cannot involve less than three coordinates (even in Minkowski space), 
instead of one for the Laplacian on Riemann spaces of constant curvature. Estimations and 
comparison techniques of stochastic analysis are the main tools we use to prove our results. 
They do not include yet a full determination of the Poisson boundary, but they suggest 

1 Note however that some physical models of diffusion in a relativistic fluid are not Lorentz invariant 
since the frame of "the fluid at rest" plays a specific role : cf [D] and its references. 
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that for general Lorentz manifolds bounded harmonic functions could be characterised by 
classes of light rays, i.e. null geodesies. 

Let us now explain more precisely the content of this article. 

We consider diffusions, namely continuous strong Markov processes. We start, in Sec- 
tion 2.1, with the flat case of Minkowski space R 1,rf , and therefore with the Brownian 
motion of its unit pseudo-sphere, integrated then to yield the only true relativistic diffu- 
sion, according to [Dl]. We get then its asymptotic behaviour, somewhat simplifying the 
point of view of [D3]. 

In Section 2.2 below, we present an extension of the preceding construction to the 
framework of general relativity, that is to say of a generic Lorentz manifold. The process 
is first defined at the level of pseudo-orthonormal frames, with Brownian noise only in 
the vertical directions, and projects into a diffusion on the pseudo-unit tangent bundle. 
The infinitesimal generator we get in Theorem 1 decomposes into the sum of the vertical 
Laplacian and of the horizontal vector field generating the geodesic flow. 

In Sections 2.3 and 2.4, we deal in detail with the Schwarzschild space, which is the 
most classical example of curved Lorentz manifold, used in physics to model the space 
outside a black hole or a spherical body. 

Using the symmetry, and introducing the energy and the angular momentum, which 
are constants of the geodesic motion, we reduce the problem to the study of a degener- 
ate three-dimensional diffusion. We then establish in Theorem 2 that almost surely the 
diffusion either hits the hole or wanders out to infinity, both events occurring with pos- 
itive probability. We prove also in Theorem 2 that almost surely, conditionally on the 
non-hitting of the hole, the relativistic diffusion goes away to infinity in some random 
asymptotic direction, asymptotically with the velocity of light. Then we prove in Theo- 
rem 3 that almost surely, conditionally on the hitting of the hole, the relativistic diffusion 
reaches the essential singularity at the center of the hole within a finite proper time, and 
we describe the limit. 

We show then in Theorem 4 that the story can be continued further : namely, the 
Schwarzschild relativistic diffusion, a priori defined till its hitting of the center of the 
hole (the essential singularity of the so-called Kruskal-Szekeres space, we also call full 
Schwarzschild space), can be extended to a diffusion which crosses this singularity. This 
extended diffusion reaches soon the restricted Schwarzschild space again, where it evolves 
as before, but maybe running backward in time. Thus such hole crossing can happen then 
again and again, but without accumulation, according to Theorem 4 below, so that the 
extended Schwarzschild relativistic diffusion is well defined for all positive proper times. 

We finally study the asymptotic behaviour of this extended Schwarzschild relativistic 
diffusion, and show mainly in Theorem 5 below that there is a unique alternative possibility 
to the escape to infinity : there is indeed a positive probability that the relativistic diffusion 
becomes endlessly confined in a spherical neighborhood of the hole, with an increasing 
velocity, and a trajectory becoming asymptotically planar, with an asymptotic shape. This 
implies the existence of an infinity of SO3 -invariant harmonic functions. 
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2 Statement of the results 



2.1 A relativistic diffusion in Minkowski space 

Let us consider an integer d > 2 and the Minkowski space R 1 '^ := {£ = £) G 
R x R d }, endowed with the Minkowski pseudo-metric (£,£) := |£°| 2 — ||£|| 2 . 

Let G denote the connected component of the identity in 0(1, d), and denote by 
M. d := {£ G R M | £° > and (£, £) = 1} the positive half of the unit pseudo-sphere. 

The opposite of the Minkowski pseudo-metric induces a Riemannian metric on M d , 
namely the hyperbolic one, so that M d is a model for the d- dimensional hyperbolic space. 
A convenient parametrization of M d is (g,0) G R+ x S^ 1 , given by g := argch(£°) and 

Q := i I y^ ! 2 — 1 . In these coordinates the hyperbolic metric writes dp 2 + sh 2 g\d9\ 2 , 

and the hyperbolic Laplacian is := 7— + (d — l)coth g — — h sh ~ 2 g x A g , A 

og 2 ag 

denoting the Laplacian of § d_1 . The associated volume measure is |sh g\ d ~ l dgd9 . 

Note that G acts isometrically on R 1,d and on M d , and that the Casimir operator on 
G induces on HP* the hyperbolic Laplacian. 

Fix a > , and denote by C a the a-relativistic Laplacian, defined on R M x M d by 

2 TUT 

that is to say C a f := (p, grad^/) + \ A^/ . This is a hypoelliptic operator. 

Given any (£o,Po) G R 1,d x M d , there exists a unique (in law) diffusion process 
(£s,p s ) , s G R+, solving the £°-martingale problem, that is to say such that for any com- 

rs 

pactly supported / G C 2 (R M xtf), /(&,p s ) - (£o,Po) - / £ a f(£t,Pt)dt is a martingale. 

J 

Note that p s is a hyperbolic Brownian motion, and that £ s = £ + / Pt dt . 

Jo 

Note also that £ s is parametrized by its arc length. Mechanically, £ s describes the 
trajectory of a relativistic particle of small mass indexed by its proper time, submitted 
to a white noise acceleration (in proper time). Its law is invariant under any Lorentz 
transformation. 

Note that if we denote by (eo, e±, .., e^) the canonical base of R 1,d , and by (e*) the 
dual base (with respect to ( , )), the matrices Ej := e <8> e* + ej <g> belong to the Lie 
algebra of G , and generate the boost transformations. Given d independent real Wiener 
processes w{ , p s = (p°,p s ) can be defined by p s := A s e , where the matrix A s G G is 
defined by the following stochastic differential equation : 



■j o dw t 
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This means in fact that the relativistic diffusion process (£ s ,Ps) is the projection of some 
diffusion process having independent increments, namely a Brownian motion with drift, 
living on the Poincare group. This group is the analogue in the present Lorentz-Minkowski 
setup of the classical group of rigid motions, and can be seen as the group of (d + 2, d + 2) 

real matrices having the form ^ with A G G , £ G R 1 '** (written as a column), 

and G (written as a row). Its Lie algebra is the set of matrices ^ with 

(3 G so(l, d) and x G IR 1 '^. The Brownian motion with drift we consider on the Poincare 

group solves the stochastic differential equation d ^ ) ~ ( q S ^ ) ° ^ ( ^ 6< Q 

where (/3 a = crX^=i Ej w i) is a Brownian motion on so(l,d). This equation is equivalent 
to dA s = A s o rf/5 s and = A s eo ds , so that (A s ) is a Brownian motion on G . On 
functions of p — Aeo , its infinitesimal generator X^=iOCe,) 2 coincides with a Casimir 
operator, and induces the hyperbolic Laplacian, so that (p s = A s e ) is a Brownian motion 
on M d , as required. 



Then it is well known that 6 S := p s j sj \p° s \ 2 — I converges almost surely in S d 1 to 
some random limit 9^ , and that p° increases to infinity. Set also g s := argch(p°) . 

The Euclidian trajectory Z(t) is defined by £ s ( t ) , where s(t) is determined by = t . 
Let us note that the Euclidian velocity — -— = s ( t ) th g s ( t ) has norm < I (I is here the 

(XL 

velocity of light). Moreover we have the following. 

Z(t) 

Remark 1 The mean Euclidian velocity converges almost surely to 6^ G S d_1 . 



Proof We have lim s(t) = +oo , so that th^ s ( t ) = Jl — {p s(t) ) 2 goes to f. Thus we 

get almost surely lim ^Jil — ) anc i ^he result follows at once, o 

t^-oo dt 

Remark 2 The scattering amplitude, id est the law of 9^ given p , is given by the 
hyperbolic harmonic measure in the unit ball of M. d (taken as model for H d ), which has 
density proportional to P(po, ■) c! ~ 1 with respect to the uniform measure of E> d ~\ P denoting 
the classical Poisson kernel of the unit ball of M. d . See for example ([E-F-LJ], case 5 = 0). 



2.2 Extension to general Lorentz manifolds 

Let us now see how the preceding construction can be naturally extended to the frame- 
work of manifolds. 

Let A4 be a (d+ l)-dimensional manifold, equipped with a pseudo-Riemannian metric 
of signature (+, —,..,—), together with an orientation and a time direction, and its Levi- 
Civita connection. 



5 



For notational convenience, T l M will always denote the positively oriented half of the 
unit tangent bundle of M. . As in the construction of Brownian motion on Riemannian 
manifolds, we have to use the frame bundle (see [M]). 

So let G(A4) be the bundle of direct pseudo-orthonormal frames, with first element 
in the positive half of the unit pseudo sphere (in the tangent space), which has its fibers 
modelled on the special Lorentz group G . Let Vj be the canonical vertical vector field 
associated with the preceding matrix Ej , and H be the first canonical horizontal vector 
field. 

Set £:=# + tEV?. 

3=1 

Let 7Ti denote the canonical projection from G(M) onto the tangent bundle T l M. , 
which to each frame associates its first element. The canonical vertical vector fields Vki 
associated with the matrices E k i := e k <g> e* — e\ ® e\ G so(l,d) , for 1 < k < I < d , 
generate an action of SO a on G(Ai), which leaves T X M. invariant and then allows the 

identification T l M = G(M)/SO d . The Casimir operator is C = £ V? - £ V k \ . 

3=1 l<k<l<d 

Note that the matrices {Ej , E ki ; l<j<d,l<k<l<d} constitute a pseudo- 
orthonormal base of so(l,d) (endowed with its Killing form). 

d 

Lemma 1 The operators H , V? , C , C do act on C 2 functions on the pseudo-unit 

3=1 

tangent bundle T X M. , inducing respectively : the vector field Cq generating the geodesic 
flow on T X M. , the so-called vertical Laplacian A v , A v again, and the generator Q : = 
Cq + ^- . 

More precisely, for any test-function F on T X M., we have on G(M) : 

(CqF) O TTi = H (F O TTx) , (A V F) O TTi = C(F O TTi) . 

Besides, in local coordinates (x\ e*) , with ej = ej : Vj = e^ -J=^ + e\ -J^ , and ((g kl ) 
denoting in these coordinates the inverse matrix of the pseudo -Riemannian metric of M.) : 

i^F) o tt! = f Vf{F o ttO = ((e fe 4 - g») ^ + de k ^)Fon l . 

Proof Let us observe that for any u = (x, e , .., e d ) £ S0 ljd M , if (u s ) denotes the 
horizontal curve such that u = u and n*u = x = e , then (7Ti(ti s )) is the geodesic 
generated by ni(u) = (x,e ) . Hence for any differentiable function F on T X A4 , we have 

H (F o 7Ti)(«) = ^F o 7Ti(u s ) = ^FfaOO) = CoFMu)) . 
ds ds 

Another way of expressing this is to recall that H commutes with the rotation vertical 
vectors V k i . It is also classical that the Casimir operator C commutes with all vertical vec- 
tors Vu , Vj . Moreover, since the rotation vectors Vki act trivially on T X M, the operators 

d 

C and V? induce the same operator on T l M. 

3=1 
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Then as eQ[e tVj (x, e)] = e§(x, e) cht + e*(x, e) slit and e k \e tY ^ (re, e)] = e^(x, e)cht + 
eg(x,e)sht, we have indeed Vj/(x,e) = jfcf[e tv i(x,e)] = e k ^f(x,e) + eg e) , id 

J 

est Vj = e^-^k + e o7Tfe • Using that ene l — J2 eke \ — 9 kl , we deduce immediately : 

d n2 d a f) d f\2 d f\2 

V\/ 2 - (e k e l -a kl ) — + V e k — - + de k — - + 2 V e k e l — + VeV — 

^ - le e j d ^ + ^ e j ^ + e ^ + ^ Q ^ + ^ e e Q ^ , 

which reduces to the formula of the statement in the particular case of a function de- 
pending only on (x, eo). In accordance with the commutation relations arguments above, 
Y,j=i ° ^i) i s a function depending only on (x, eo), id est a function on T X J\A. o 

Now, according to Section 2.1, the relativistic motion we will consider lives on T X J\A 

2 

and admits as infinitesimal generator the operator Q = C + ^ A v of Lemma 1 above. 

If M. is the Minkowski flat space of special relativity, it coincides with the diffusion 
defined in Section 2.1 above. 

To construct this general relativistic diffusion, we use a kind of stochastic development 
to produce a stochastic flow on the bundle G(J\4), as is classically done to construct the 
Brownian motion on a Riemannian manifold. But we have now to project on T l M. and 
no longer on the base manifold J\4, and to put the white noises on the acceleration, id est 
on the vertical vectors, and no longer on the velocity, id est on the horizontal vectors. 

To proceed, let us simply fix \l/o £ G(J\4) and a M. d -valued Brownian motion w = (w J s ) , 
and let us consider the G(J\4) -valued diffusion ^ = (\l/ s ) G G(M) solving the following 
Stratonovitch stochastic differential equation : 

(*) ^ s = + f H Q {%) dt + a [ S J2 Vj(%) o dw{ . 
Jo Jo ^— [ 

By Lemma 1, the stochastic flow defined by (*) commutes with the action of SOd on 
G(J\4), and therefore the projection (£ s ,£ s ) := (£ s , e o( s )) — 7r i(^'s) defines a diffusion on 
T l M. ; namely this is the relativistic diffusion we intended to define and construct. 

The following theorem defines the relativistic diffusion (£ s ,£ s ), possibly till some explo- 
sion time. The vector field C denotes the generator of the geodesic flow, which operates 
on the position ^-component, and denotes the vertical Laplacian (restriction to T l M. 
of the Casimir operator on G(j\4)), which operates on the velocity ^-component. 

Theorem 1 1) The G(J\4) -valued Stratonovitch stochastic differential equation 

d 

(*) d^ s = H (V S ) ds + a J2 VjiVs) o dw{ 

3=1 

defines a diffusion (£ s , £ s ) := 7Ti(^ s ) on T l J\A, whose infinitesimal generator is C + ^- A v . 
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2) If £ (s) : T^ s M. — > T^ M. denotes the inverse parallel transport along the C l curve 
{Cs' I < s' < s), then ( s := £ (s) £ s is an hyperbolic Brownian motion on T^ M. . 

Therefore the path (£ s ) is the development of a relativistic diffusion path in the Minkowski 
space T^ Q M. . 

Remark 3 In local coordinates (x l ,e*) , with ej = e k -£^ , = (£s ; e (s), ..,ed(s)^j , 
T k t denoting as usual the Christoffel coefficients, the equation (*) writes : 

d 

dx\ = ei(s)ds ; de k (s) = -r£(&)ej(s)da;* + l {jm cre k (s) o dw{ + l {j=0 }^E e i( s ) ° dw l » 

i=i 

or equivalently in the Ito form : 

d 

dx k = e k (s) ds ; de k (s) = ~T k l ^ s )e l (s)dxl + aJ2e k (s)dwl + ^e k (s)ds, and 

i=i 

dej(s) = -r£(&) ej(s) dx* + a e k (s) dw ] s + ^e k (s)ds for j > 1 , < k < d . 
Note that the martingales in the above equations for e (s), that is to say the differentials 

d 

dM k := o-^2 e k (s)dw J s , < k < d , have the following quadratic covariation matrix : 

j'=i 

K M. = (_dM^) = a 2^ is y M = ^ (e k (s)e l (s) - g k %)) , 

id est K s = (T 2 (e (s) 'eo(s) — <7 -1 (£ s )) , in accordance with Lemma 1. (Here 'e denotes 
the transpose of the column-vector eo , and g~ l denotes the inverse matrix of the pseudo- 
metric.) K s has rank d : we have indeed Kg eo = (eo *eo — Q 1 ) 9 e o — , since *eo eo = 1. 

Note that does not depend on the other frame vectors ej (j > 1), proving again 
that the projection (£ s ,£ s ) is a diffusion on the tangent bundle, as Theorem (1,1) asserts. 

Proof of Theorem 1 1) does not need any further proof. 

2) The process (( s = £ (s) £ s )s>o is continuous and lives on the fixed unit pseudo- 
sphere T£ M ■ Let £ (s) : Tg .M — > T ?s .M denote the parallel transport along the C 1 
curve (£ s / | < s' < s), and recall that 

£ 7(*)5 = - T(s) k TL(6) C , which implies that £ T(*)5 = 0,(6) C 5 

indeed 

(£ T(«)5) T(«)5 = TW5 ?(*)? rL&) c = 7(*)J Tw; rLfc) c , 

whence 

£ t(s)i = t(s)l T(s) k t(s% rijQ 4 m = t(s% rj m (c s ) 4 m . 
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Recall then from Remark 3 that (for < t < d ) 

di e s = o 4(s) o dw k s - r e jk (t s ) ii g ds . 
k=i 

Therefore we get 

d£ = (orfT(^)4 £ + T(^o^f 

= t(sT q TLfc) C ds + a tisfe £ e{{s) o ^ - <f ^(6) g 4" * 

fc=i 

fc=l fc=l 
where efc(s) := £ (s) efc(s) , for 1 < k < d and s > . 

Similarly, we have de e k (s) = a £ (s)^ £f o rfwf 1 , id est dek(s) = a( s o <iu^ . 

Observe that, for any s > , ei(s), .., ed(s)) constitutes a pseudo-orthonormal basis 
of the fixed tangent space T^ Ai . Hence, owing to Section 2.1 and Remark 3, we find 
that the velocity process (. defines a hyperbolic Brownian motion on the hyperbolic space 
T^M, isometric to M. d . 

In the reverse direction, we have of course £ s = £ (s) ( s , meaning indeed that we 
recover the C 1 curve £. as the deterministic development of the flat relativistic diffusion 
JoCds. o 

Remark 4 The equation (*) can be expressed intrinsically, in Stratonovitch or in Ito 
form, by using the covariant differential D , which is defined in local coordinates (x l , e k ) 
by (Dej) k := de k + Y k u ej dx % , for < j, k < d . The equation (*) is indeed equivalent to : 

d 

— e o( s ) ; De (s) = cr^e^s) o dw{ ; Dej(s) = a e (s) o dw{ for 1 < j < d . 

3=1 

See for example ([B], page 30) or ([Em], page 427). 

2.3 The restricted Schwarzschild space So 

This space is commonly used in physics to model the complement of a spherical body, 
star or black hole; see for example [DF-C], [F-N], [L-L], [M-T-W], [S]. 

We take M = S := = (f, r,d)eRx [R, +oo[x§ 2 } , where R G R+ is a parameter 
of the central body, endowed with the radial pseudo-metric : 

(1 - f ) di 2 - (1 - f )~ 1 ^r 2 - r 2 |d(9| 2 . 

The coordinate t represents the absolute time, and r the distance from the origin. 
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In spherical coordinates 9 = (</?, ip) G [0, 7r] x (R/27rZ), we have (c^l 2 = dip 2 + sin 2 ^^ 2 . 
The geodesies are associated with the Lagrangian , where 

2 L(£ = (1 - ?) i 2 " (1 " f )^ r 2 - r V - r 2 sin 2 v ^ , 
and the non- vanishing Christoffel symbols are : 

^ = -Kr = ; n t = ; r w = a-r ; = (R-r)sm*<p ; 

The Ricci tensor vanishes, the space <5>o being empty. A theorem of Birkhoff (see [M-T-W]) 
asserts that there is no other radial pseudo-metric in Sq which satisfies this constraint. 
The limiting case R = is the flat case of special relativity, considered in section 2.1. 



2.3.1 The stochastic differential system in spherical coordinates 

Let us take as local coordinates the global spherical coordinates : 

z=(e,e,e,e) -.= (t,r,^). 

According to Remark 3, the system of ltd stochastic differential equations governing 
the relativistic diffusion (£ s ,£ a ) writes here as follows: 

dt s = e (s) ds , dr s = ej(s) ds , d<p s = el(s) ds , dip s = el(s) ds , 

del{s) = 3 4 e° (s) ds - e° (s) ej( a ) ds + dM° , 

del{s) = "4 el(s) ds + e£(s) 2 ds - e° (s) 2 ds + (r s - R) e 2 (s) 2 ds 

+ (r s - R) sin 2 <p a e 3 (s) 2 ds + dM) , 

del(s) = 4 e o( s ) ds - j- ej(s) e\{s) ds + sin (p s cos y2 s eo(s) 2 rfs + dM 2 , 
rfeo(s) = 4 e o( s ) ds - £ej(s) ejj(s) ds - 2cotg</? s e 2 ,(s) ejj(s) ds + dM 3 , 

where the martingale M s := (M°, Mj, M 2 , M^) has the following rank 3 quadratic covari- 
ation matrix : K s = a 2 (e (s) 'eo(s) — <7 -1 (£ s )) • 



2.3.2 Energy and angular momentum 

We shall use widely the angular momentum b := r 2 9 A 9 , 
the energy a :— (1 — j) i — ^ , and the norm of 6 : 6 := |6| = r 2 U , with [/ := \9\ . 
Set also T :— f , and accordingly 

T s := r s = ej(s) , U s := \9 S \ = \J el(s) 2 + sin 2 ip s efj(s) 2 , and D := min{s > | r s = R} . 

Standard computations yield the following : 
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Proposition 1 1) The unit pseudo-norm relation (which expresses that the parameter 
s is precisely the arc length, id est the so-called proper time) writes 

T 2 = a 2 s -(l-R/r s )(l + b 2 Jr 2 ). 

2) The process (r s , o s , b s , T^) is a degenerate diffusion, with lifetime D, which solves 
the following system of stochastic differential equations : 



b 2 R 

dr s = T s ds , dT s = dMj + ^T s ds + (r s - \R) -±ds - — ds , 



s 



K := a 2 



da s = dM a s +^a s ds , db s = dM b s + ^-b s ds + ^— ^ ds , 

z o s 

with quadratic covariation matrix of the local martingale (M a , M b , M T ) given by 

i 2 s - 1 + £ a s b s a s T s 

a s b s b 2 s +r 2 s b s T s 
\ a s T s b s T s T 2 + l-f. 

We get in particular the following statement, in which the dimension is reduced. 

Corollary 1 The process (r s ,b s ,T s ) is a diffusion, with lifetime D and infinitesimal 
generator 

8 a 2 d 2 a 2 8 d 2 



dr 2 * ' db 2 2b * ' db db&T 

<j 2 / „ „> d 2 /3<r 2 , . 1 ! R > d 



We have the following result on the behaviour of coordinate a s . 

Lemma 2 There exist a standard real Brownian motion w s , and a real process i] s , 
almost surely converging in R as s S D , such that a s = exp(<r 2 s + aw s + i] s ) for all 
s E [0,D[. In particular a s almost surely cannot vanish, which means that time t s is 
always strictly increasing. 

Proof Proposition 1 above shows that (a 2 — 1) a 2 ds < (dM®) < a 2 a 2 ds , for < s < D . 
So that we have almost surely (as s — > oo , when D = oo) : 

f S f S f s 

loga s -loga = 3(j 2 s/2-i / a7 2 (dM?) + / a' 1 dM? > a 2 s+ aT 1 dM^ 

Jo Jo Jo 

= a 2 s + oi y J\- 2 {dM^)) = a 2 s + o{s) . 

r-D 

Since (1 — y-) < a 2 , this implies / (1 — ^) a J 2 ds < oo almost surely. 

J 
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Consider then the standard real Brownian motion w defined by 



dMg = a y a? s — (1 — R/r s ) dw s , and the process 77 defined by the formula in the statement. 
We have : 



drj s = d(\oga s ) — a 2 ds — a dw s = \ (1 — R/r s ) a s 2 a 2 ds + (^J 1 — (1 — R/r s )a~ 2 — l) a dw s , 
and then for any s < D : 



.^o^Al- R/n) a' 2 dt-af pm « 
2 Jo Jo 1 + A _d_ 



2 



+ Jl-(l-R/r t )a 



-2 



dw t , 



which converges almost surely to a finite limit as s /* D , since almost surely for all 
s E]0,D[ : 

2.3.3 Asymptotic behaviour of the relativistic diffusion (^ s ,^ s ) 

We see in the appendix (Section 4.1) that in the geodesic case a = five types of 
behaviour can occur, owing to the trajectory of (r s ) ; it can be : 

- running from R to +oo, or in the opposite direction ; 

- running from R to R in finite proper time ; - running from +oo to +oo ; 

- running from R to some R\ or from R\ to +oo, or idem in the opposite direction ; 

- running for endlessly in a bounded region away from R. 

More detailed results can be found in [L-L] and mainly [M-T-W]. A full treatment is 
given, for future reference, in the appendix below (Section 4). 

The stochastic case a ^ can be seen as a perturbation of the geodesic case a = ; 
however the asymptotic behaviour classification is quite different. 

Theorem 2 1) For any initial condition, the radial process (r s ) almost surely reaches 
R within a finite proper time D or goes to +oo as s — > +oo (equivalently : as t(s) — > +oo 
if Oq > 0, and as t(s) — > — oo if Oq < 0). 

2) Both events in 1) above occur with positive probability, from any initial condition. 

3) Conditionally on the event {D = oo} of non-reaching the central body, the Schwarzschild 
relativistic diffusion (£ s ,£ s ) goes almost surely to infinity in some random asymptotic di- 
rection of M 3 , asymptotically with the velocity of light. 

Note in particular that the relativistic diffusion almost surely cannot explode before 
the finite proper time D . 

The proof we give for this theorem is rather long. It is postponed till Section 3 below. 
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2.4 The full Schwarzschild space S 

The full Schwarzschild space S, also known as the Kruskal-Szekeres space (see [DF-C] , 
[F-N], [S], and especially [M-T-W] and its historical account page 822), can be defined by 
extending the previous restricted Schwarzschild space S as follows. On S , set 



u 



:= x e r/{2R) x ch (^) and v := x e r/{2R) x sh (^) . 



Note that (— — l) x e r ^ R = u 2 — v 2 , and that the Schwarzschild pseudo-metric ex- 
presses in the Kruskal-Szekeres coordinates (u,v,9) as: 

(**) i5! e r ' R (dv 2 -du 2 )-r 2 \dB\ 2 , 
r 

where r = r(u 2 — v 2 ) , r denoting the inverse function of [r i— > (-^ — 1) e r / fl ] (which is an 
increasing diffeomorphism from R + onto [— l,+oo[). 

In those Kruskal-Szekeres coordinates, we have S = v,6) G I 2 x § 2 u > \v\j. 
The full Schwarzschild space S is now defined as 

S ■= {(u, v, 6) G R 2 x § 2 | u 2 - v 2 > -l} , 

and is equipped with the pseudo-metric defined by (**) above and by r = r(u 2 — v 2 ) . 

S contains So , — So (isometric to «So), two isometric copies of the hole: 

Tt : = {(«, v , 6) G S v > and — H , 
and the boundary between ±So and ±7i , which is {r = i?} = {u = ±v }. 

It is a Lorentz manifold, to which our general construction of section 2.2 applies. 

The energy and angular momentum are extended to T X S by setting 

2R 2 



a := 



r 



e r/R (uv - vu) and b := r 2 6 A 6 . 



As before we set T := f = 2 ^- e r / R (uii — vv) and b := |6| = r 2 U . Recall that 
a and b are constant along geodesies. The unit pseudo-norm relation writes as before : 

« 2 -r 2 + (f-D(^ + i)=o, 

T € r ^ R / b 2 \ 

or equivalently u 2 — v 2 H — ( — + 1 ) = , which implies \v\ > \u\ , whence v > \u\ 

along a timelike path. The following correspondences between a line-element (£, £) G T 1 ^ 
and its projection ^ G 5 are easily deduced : (£, £) G {r > R ; a > 0} £ G iS>o , 

E {r > R; a <0}&£ e -So, (£, G {r < R ; T < 0} & £ G -W , 
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(Z,t)e{r<R;T>0}&SeH, (f,0 G {r = R; a = T = 0} & £ G = v = 0}. 



Now two other coordinates, namely the so-called inward and outward Eddington- 
Finkelstein coordinates u~ and u + , prove to be very convenient for performing calcu- 
lations on S . They are defined, not on the whole S , but u~ on S fl {u + v ^ 0} and 
u + on 5 fl {ii - w ^ 0}, by : 

u~ := 2R log |w + f | and u + := — 2R log |u — t> | . 
In those Eddington-Finkelstein coordinates, the metric expresses as : 

(l-^){du ± f±2du ± dr -r 2 \d6\ 2 , 
and the energy expresses as 

a = (1 - 5)u± ±r = . 

We shall need (from Section 2.4.2) to complete the space S, in S := S U dS, where 

OS := { ± (u,v,6) | (u,u,0) G M 2 x § 2 ; m 2 - v 2 = -l} , 

meaning that we identify the opposite points above the singularity {r = 0}. Note that 
r, u + , u~ are naturally continued to S, that r = ^=^> v 2 = u 2 + 1 <^=^> u + = u~ in S , 
and that 

w+ + r + i? log |^ - 1| = u~ - r - R\og\^-l\ on 5 n {\u\ ^ \v\} , 

this quantity being equal to t on ± So . The R- valued absolute time t is continuously 
extended to S fl {u = v = 0} c , by setting t = +oo on {u = v ^ 0}, t = — oo on 
{u = —v 7^ 0}, t = 2Raxgth(u/v) on ±7Y, and t = 2Raxgth(v/u) on — S$ (as on «So). 

Note that the region A := {u = v = 0} appears right away as exceptional, as the only 
part of S where t cannot be continued, and the only part of S where both u^ 1 explode. 



In the Kruskal-Szekeres coordinates, a path (u s ,v s ,9 s ) is timelike if and only if v s > 
| . This implies that any timelike path started in H has to hit {v = Vu 2 + 1 } C {r = 0}, 
and that any timelike path started in —H has to hit {r = R}. In particular, a timelike 
path started in — Ti and avoiding A has to enter the region {r > R} through {t = — oo}, 
appearing either as a particle born and then evolving in So (for ever, or entering then Ti 
through {t = +oo}), or the analogue through — So , which could be viewed as the case of 
an antiparticle through S . 

These dynamics of the full Schwarzschild space S show that the inward Eddington- 
Finkelstein coordinate u~ is appropriate to the study of timelike paths started in Sq U 7i, 
till they hit {r = 0} (and even to extend them further, see section 2.4.2 below), and that 
the outward Eddington-Finkelstein coordinate u + is appropriate to the study of timelike 
paths started in — Ti and entering ±S , till they hit Ti. 



14 



2.4.1 Diffusion in the full space S : hitting the singularity 

We follow the same route as in the restricted Schwarzschild space, to express the rel- 
ativists diffusion on T 1 ^. Let us proceed, using the Eddington-Finkelstein coordinates 
u ± . The Lagrangian writes 

2 L(£, = (1 - f) {i^f ± 2 r - r 2 2 - r 2 sin 2 y? V> 2 . 

We apply then Remark 3, to get the Ito stochastic differential equations of the rela- 
tivistic diffusion in the full Schwarzschild space : 

diif = a dMf + ^f-iifds±£s (iif) 2 ds =p r s (0 2 + sin 2 Lp s ip 2 s ) ds ; 
dr s = a dM r s +^r s ds+{^-l)A{uf) 2 dsT$ufr s ds + {r s - R) (y? 2 + sin 2 Vs if 2 s ) ds ; 

£t I S Zl ' S 1 S 

d(p s = a dM^ + ^ <p s ds — -jr f s ip s ds + sin ip s cos (p s %l) 2 ds ; 
dip s = a dMf + ^ ip s ds - £ r s ip s ds - 2 cotg (p s (p s ip s ds , 

for some continuous local martingale (M ± , M r , Mf , Mf) , having quadratic covariation 
matrix (according to Lemma 1 and Section 2.3.1) : 

/ (iif) 2 iif r s =f 1 uf <p s iif f s \ 

iif r s =F 1 r 2 s + 1 - f- B f s ip s r s ip s 

uf (p s r s (p s (f 2 + r- 2 _ Lf> s ip s 

V ufip s r s ip s (p s ij s ij 2 + (r s sin ip s y 2 / 

We have again a reduced diffusion (r s ,b s ,T s ) (with minimal dimension), solving the 
same system of Ito stochastic differential equations as before. 

This system of stochastic differential equations has been derived using Eddington- 
Finkelstein coordinates, so that it is valid a priori outside {u = v = 0}. But the smooth 
functions (r, a, b, T) of the relativistic diffusion have an Ito decomposition with continuous 
coefficients, so that the formulas involving them hold without restriction. 

From the pseudo-norm relation, we see that T s cannot vanish in the region {r < R}. 
As r s enters this region necessarily almost surely with derivative < (indeed \Tu\ = 
\cle>\ > by Lemma 2), r s is then necessarily strictly decreasing. Precisely, we have the 
following. 

Theorem 3 The relativistic diffusion in T X S either escapes to infinity, or enters above 
Ti at time D and converges to the singularity within some finite proper time D' . More- 
over, in the second case we have almost surely : 

1) for s > D, r s decreases and hits at proper time D' , with D < D' < D + | R ; 
moreover lim S/ /<£>/ T s = — oo ; 

2) 9 S , uf , b s and a s converge to finite limits as s /* D' , and b D i cannot vanish; 
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3) as s y D' , we have the following equivalents: 

r a ~[\b D ,\fR{D' - s)] 2/5 and T s ~ -b D , VR x [| b D ,^R (£>' - s)]~ 3 / 5 . 

Remark 5 The equivalents in Theorem (3,3) above can be specified further. Indeed, 
we have almost surely, as s \ : 

rjy-. = [f b D ,VRsf 5 x (l + S ) 2 / 5 + 0(^log|log S |)) , 

and 

T^_ s = -b D ,VR x [| fo^v^s]- 3 / 5 x (l - 2(§£ s) 2 / 5 + 0(0* log | log s|)) . 

Indeed, using the stochastic differential equation of T s and the iterated logarithm law, 
together with the equivalents in Theorem (3,3), we deduce easily these more precise asymp- 
totic expansions near D' . 

Remark 6 We know from Theorem 1 that the relativistic diffusion can start from any 
initial condition in the full space T l S . When it starts above — H , the pseudo-norm relation 
forbids any vanishing of T s (which has then to remain > ), till the level {r = R} is hit, 
which takes a proper time less than -rrR/2 , for the very same reason as in Theorem (3,1). 
When the diffusion starts above {r = R}, it enters {r ^ R} at once, as any timelike path. 
Note that above A = {u = v = 0} C {r = R}, we have necessarily T = a = . Moreover 
it can be proved that T X A is polar for the relativistic diffusion. So, when starting above 
—H, the relativistic diffusion enters then above ±<So , before possibly entering later above 

n. 

We postpone the proof of Theorem 3 till Section 3 below. 

A part of this proof is based on the following proposition, which allows to recover the 
whole relativistic diffusion from the reduced relativistic diffusion (r s ,b s ,T s ). 

Proposition 2 The spherical coordinate 6 S satisfies the following stochastic differential 
equation (conditionally on the reduced relativistic diffusion (r s , b s , T s )) : 

for some standard real Brownian motion (3. , which is independent of (r.,b.,T.) . 
Moreover, 9 S /U S converges in § 2 as s / D' , almost surely. 

We postpone the proof of this proposition till Section 3 below. 

Corollary 2 The curve in the full space S defined by the image of the trajectory 
{(r s , uj, 6 S ) | s < D'} admits almost surely a semi-tangent at the center of the hole. 
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Proof Using the strict monotonicity of r s near the singularity, we see that it is sufficient 
to verify the left-differentiability of the curve (r s h- > (u~,r s 9 s )) at s = D'. Now using 
Theorem (2,3), as s/D',on one hand we have 

du~ /dr s = u~ /T s = (a s /T s + f ) r s / (r s — R) ~ —r s /R — > , and on the other hand 

d{r s 9 s )/dr s = 9 S + {9 s r s /T s ) -> D , G § 2 , since 

\9 s r s /T s \ = b s /(r s T s ) ~ -iT 1 ^ [| b D ,^R{D' - s)} 1 / 5 -> . o 

Remark 7 Since we have |# s | = C/ s = &s/r^ , we observe the explosion of the spherical 
speed # s , as well as of the radial speed T s , and as of the speed iij, at proper time D', id est 
at the hitting of the singularity {r = 0}. However we also just saw that (r s , b s , a s , 9 S , u~) 
is left-continuous at D' , and that moreover the curve in the space S defined by the image of 
the trajectory {(r s ,uj ,9 S ) \ s < D'} admits almost surely a tangent at any point. Indeed 
it happens that the explosion of the derivatives does not forbid to define a continuation 
of the relativistic diffusion after the finite hitting proper time D'. This will be indeed the 
purpose of Section 2.4.2 below. 

2.4.2 Regeneration through the singularity: entrance law after D' 

We see from Theorem 3 and Remark 7 that the set of endpoints of the relativistic paths 
(at proper hitting time D' of {r = 0}) identifies with the boundary dS, which we defined 
in Section 2.4, identifying pointwise the outward &H and inward d(—7i) boundaries. Thus 
we have indeed dS = d^S. 

In this identification, a differentiable inward path ending at {r = 0} can be continued 
by a differentiable outward path, so that the M 3 -valued curve r9 is differentiable at any 
point. In particular, geodesies are thus well defined for any proper time, and there are 
geodesies which cross endlessly the singularity {r = 0}, namely those which are described 
in case 1.2 (and are met also in cases 2.2.1, 2.2.2, 2.5.2, and 2.6), completed by Remark 13, 
of Section 4. For generic values of parameters (a,b,k), such geodesies are dense in some 
disk of R 3 (centred at 0). 

Since the diffusion can hit the singularity {r = 0} in finite time, it is natural to look 
for an entrance law, allowing to continue it after time D' . Clearly it as to enter above 
—7i. Thus we have to define for the diffusion on T 1 ^ a family of entrance laws above the 
singularity {r = 0}, and more precisely on the boundary dS. 

Let Q denote the generator of the relativistic diffusion, acting on C 2 {T X S) . Theorem 
3 allows us to extend the relativistic diffusion to a continuous strong Markov process on 
T X S U dT x S, provided we establish the following proposition. 

Proposition 3 The martingale problem associated with Q has a unique continuous so- 
lution, starting from any point of dS. 

We postpone the proof of this proposition till Section 3 below. 
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2.4.3 The relativistic diffusion for all positive proper times 

The existence and uniqueness of the entrance law (in Proposition 3) allows to prove the 
first assertion of the following. 

Theorem 4 There exists a unique continuous strong Markov process on S with positive 
lifetime inducing the relativistic diffusion on S . 

The lifetime of this extended relativistic diffusion is almost surely infinite. 

For such a process, we define an increasing sequence of hitting proper times D k as 
follows. Let Dq := D G [0, oo] denote the hitting time of {r < R}, Di := D' denote the 
hitting time of {r = 0}, and set by induction, for any n G N* : 

D 3n := inf{s > D 3n -i \ r s = R} , D 3n+1 := inf{s > D 3n \ r s = 0} , and 

D 3n+2 := inf{s > D 3n+1 \ r s = R}. Finally consider := sup n D n . 

This is obviously an increasing sequence of stopping times, strictly increasing as long as 
it is finite. The preceding section 2.4.2 extends in a unique way the law of the relativistic 
diffusion to the proper time interval [0, [. It is clear that the process cannot be extended 
continuously beyond . 

This proves the first assertion of Theorem 4, within the lifetime . We postpone till 
Section 3 the proof of the second assertion of this theorem : Doq is almost surely infinite. 

Remark 8 We saw in Theorem (3,1) that D\ < D + \ R- Now exactly the same 
reason shows that D 3n+ \ — D 3n and D 3n+ 2 — D 3n+ \ are < 7r R/2 , as long as these times 
are finite. The time intervals [D 3n _i, D 3n ] correspond to the excursions outside the hole, 
and the time intervals [D 3n , D 3n+2 \ correspond to the excursions inside the hole. Moreover 
we see that D k becomes infinite if and only if k = 3n and the process escapes to infinity 
during its n-th excursion outside the hole. 

Recall that (according to Lemma 2) during every excursion outside the hole {r < R}, 
the diffusion can have its absolute time coordinate t — t s strictly increasing from — oo to 
+oo or strictly decreasing from +oo to — oo (this case can be seen as the antiparticle case), 
depending on the sign of a at the exit of the hole. 

Moreover, the M 3 -valued curve (s i— > r s 9 S ) is differentiable at any proper time s, 
whereas the M- valued curves (s i— > uf) present a cusp (with half a tangent : this appears 
in the proof of Theorem 3, where we saw that ii ~ —rT/R — > ±oo near D') at proper 
times D 3n+ i (and are differentiable at any time s ^ D 3n+1 ). 

Remark 9 The Liouville measure is invariant for the extended relativistic diffusion on 
S. It induces the invariant measure drdadT of the autonomous diffusion (r s ,a s ,T s ). 
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2.4.4 Capture of the diffusion by a neighbourhood of the hole 

The preceding section leads naturally to the following question : can the extended 
relativistic diffusion cross infinitely many times the hole, as some geodesies do ? 

Note first that it was clear from the preceding sections 2.3.3, 2.4.2, and 2.4.3, that 
there is, for any n G N and any initial condition, a positive probability that the extended 
relativistic diffusion crosses exactly n times the hole and thereafter goes away to infinity. 
So the following theorem asserts essentially that the limiting case of n = oo crossings of the 
hole can also happen, thereby completing the picture of all possible asymptotic behaviours 
of the extended relativistic diffusion. 

Moreover it says that this last case corresponds to an asymptotic confinement of the 
relativistic diffusion in the vicinity of the hole. 

Theorem 5 Almost surely, from any initial condition, the extended relativistic diffusion 
can have only two types of asymptotic behaviour, each occurring with positive probability : 
Either 

1) r. and \a.\ go away to infinity, and 9. converges; 

or 

2) g := limsup r s G [R, 3R/2] , liminf r s = , b. goes away to infinity, 



Moreover 

(i) For any e > 0, if r > 3R/2 and if T is large enough, then the probability that 
the relativistic diffusion goes away to infinity is at least 1 — e. 

(ii) For any e > 0, if r <E]R,3R/2[ , T = , and b is large enough, then 



The proof of this theorem is rather delicate. It will be presented in Section 3 below. 

Remark 10 For any open interval ]x,y[c]R,3R/2[, F.(x < g < y) is a non-trivial 
SO3 -invariant (£ + ^-A^-harmonic function on T 1 M . 

The following result describes more precisely what happens when the diffusion is cap- 
tured by a neighbourhood of the hole : while (according to Theorem 5) they are asymp- 
totically planar, they exhibit progressively another type of regularity. 

Corollary 3 In the second case of Theorem 5, and more precisely conditionally on the 
event g < 3R/2, the times D' n of the first maxima of the radius r. at each excursion out 




F(\g-r \ <£)>!-£. 




converges as n 



00, towards 




dr 



(which is well defined if and only if g < 3R/2). 



s J[R-r + £ 2 r 3 ] r 
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The proof of this corollary will close Section 3 below. 

Remark 11 The result of Corollary 3 concerns the time intervals [D^ n+ i, D' n ], that is to 
say the upcrossings from the singularity to the successive tops of the limiting trajectories. 
It is very likely that the same result is valid as well for the downcrossings, that is to say 
the time intervals [D' n _ 1 , D 3n+1 ], yielding the same angular random limit (the sign of T s 
compensating for the interchange of the bounds D' n _ 1 and D 3n+1 in the integral). 

So another statement in the spirit of Corollary 3, but which demands some more work, 

fD'n 

should be : almost surely lim / d6 s = ± 2 \1> , where 

_ re dr _ r 1 dr 

Jo yV [i?-r + £ 2 r 3 ] Jo [(1 - r 3 )(R/g) - (1 - r 2 )r] 

is a strictly increasing continuous function of g/R, from [1, §[ onto [|, oo[ . 

It is thus likely that the shape of the excursions should approach more and more the 
null geodesies, id est the light rays. See the appendix, Section 4.2. 



3 Proofs 

3.1 Proof of Theorem 2 

In the proof of this theorem, we shall use the following very simple lemma. 

Lemma 3 Let M. be a continuous local martingale, and A. a process such that 
liminf A S /(M) S > almost surely on {(M)^ = oo}. Then lim (M s + A s ) = +oo 
almost surely on {(M)^ = oo}. 

Proof Writing M s = W((M) S ) , for some real Brownian motion W, we find almost 
surely some e > and some sq > such that A s > 2e(M) s and \M S \ < e(M) s for 
s > s . Whence M s + A s > e (M), for s > s . o 

We prove now successively the 3 assertions of Theorem 2. 
1) Almost sure convergence on {D = oo} of r s to oo. 
This proof will be split into six parts. 

Let us denote by A the set of paths with infinite lifetime D such that the radius r s 
does not go to infinity. We have to show that it is negligible for any initial condition 
x = (r, b, T) = (r , bo, T ) belonging to the state space [R, oo[xR + x R . 

The cylinder {r = 3i?/2} plays a remarquable role in Schwarzschild geometry. In 
particular, it contains light lines. We see in the following first part of proof that we have 
to deal with this cylinder. 
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(i) r s must converge to 3R/2 , almost surely on A. 

Observe from the unit pseudo-norm relation (Proposition 1, 1) that \T s \/a s is bounded 
by I. Let us apply Ito's formula to Y s :— (1 — |^-) T s /a s : 

with some local martingale M having quadratic variation : 

(dM s) = (i-m^-f-X 1 -^^ < - 2 -s 2 ds. 



Now |Yg| is also bounded by 1. Hence Lemma 2 implies that the last two terms in the 
expression of Y s above have almost surely finite limits as s — > oo . Idem for (M s ) , and 
then for M s . Moreover the two remaining bounded variation terms in the expression of Y s 

above increase. As a consequence, we get that Y s , / — ^ dt , and / (1 — f^) 2 — ^ dt 

Jo a t rf Jo 2rt a t r$ 

dt 



a t r 2 ' 



converge almost surely in R as s — > oo . So does also / 

Jo 

■ i a , R , b 2 In _ 2 T 2 b 2 1 , , 

Now using that — < aH — - H — r = — - H H by the unit pseudo- 

r i \ r \ ar 2 a j j ar 2 ar 4 ar 2 

norm relation, we deduce that almost surely 

f < X ~ S» ! IS <Vr.)| * = fd " » 2 f * S jf (1 - g)' | * <» . 

This implies the almost sure convergence of (l ~ |~ + fpr) / r s > an d therefore of (l/r s ) . 
Since lim (l/r s ) cannot be on A , we have necessarily lim r s = 3R/2 almost surely 

on A , from the convergence of / (1 — f^) 2 -77 dt . 

Jo v 2rt> rj 

(ii) b s /a s converges to 3R\/3/2, and T s /b s goes to 0, almost surely on A. 
Indeed, Ito's formula gives (for some real Brownian motion w) 



r 2 1 — — /■« r 2 r s R b 2 

rs dw s + 2a 2 ^ds- 3a 2 / (1 - -) % ds . 



a 2 Oq Jo a 2 s \ b 2 s a 2 s Jo a 2 Jo r s a A s 

Since by the unit pseudo-norm relation we have ^ < t 2 /(1 — , whence b 2 /a 2 bounded 
on A, the above formula and Lemma 2 imply the almost sure convergence of b 2 /a 2 on 
A . Indeed the bounded variation terms converge, and as b 2 /a 2 is positive, the martingale 
part has to converge also. Using the unit pseudo-norm relation again, we deduce that 
rp2 b 2 1 

— |- = 1 — (1 — H — o ) has also to converge, necessarily to 0, since otherwise we 

aj Ts K ajrj a%> 
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would have an infinite limit for T s , which is clearly impossible on A. The value of the limit 
of b s /a s follows now directly from this and from (i). 

pco />oo 

(iii) We have almost surely on A : / (r t — ^r) 2 b\ dt < oo , and / T 2 dt < oo . 

JO JO 

Let us write Ito's formula for Z s := (r s - ^) T s = \ ± (r s - ^) 2 : 

Z s = Z + M s + (r s - r )(r s + r - 3R) + [' T 2 dt + f\r t - f ) 2 b 2 % - f [\r t - f ) % , 

Jo jo rf Jo ri 

where M. is a local martingale having quadratic variation given by : 

(M) s = a 2 J\ rt -f) 2 (l-^ + T 2 )dt. 

rs . 

Note that if (M)^ = oo , then by (ii) above hm / (r t - ^) 2 b 2 — j (M)< 

Jo 



oo . 



Note moreover that in this case / \r t — ^rl < \ [ (r t — ^r) 2 b 2 —rX\[ 777 is also 

Jo 2 r 2 _ V Jo r t \ Jo b 2 t 

negligible with respect to / (r t — ^) 2 b 2 

J 



n 

On the other hand, we must have liminf \Z S \ = on A . 

s^oo 

Therefore we deduce from Lemma 3 that necessarily (M)^ < 00 , and then that M s 
has to converge, almost surely on A . 

Using again that liminf \Z S \ =0, we deduce the almost sure boundedness and con- 



s^oo 

/•oo r-oo d£ 

vergence on A of J T 2 dt and of J (r t — ^) 2 b 2 — . 



(iv) (r s — 4^) b s and T 2 /b s go to as s — > 00 , almost surely on A . 

Indeed, on one hand we deduce from (iii) that (for some real Brownian motion W.) 

(r s -f) 2 b s = aW[n rt -fm + r 2 )dt]+2[\ rt -f)TAdt^ 

J J « 

has to converge almost surely on A as s — > 00 , necessarily to since it is integrable with 
respect to s . 

On the other hand we have for some real Brownian motion W', by It 6 formula : 



rp2 rp2 pS rpA rp4 rp2 rs rp2 rS Ji 

Jo rfbt Jo rt bt 1 Jo bf 



Recall from (i) that ^ — > and that b t ~ 3R ^ a t . Thus using (iii) we see easily that all 
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integrals in the above formula converge. Hence we deduce the almost sure convergence of 
s i— > T 2 /b s on A , necessarily to , since it is integrable. 



(y) It is sufficient to show that J \r t — 
almost surely on A . 



J \r t -™ \ \T t \ b\ dt < oo , and that J T? dt < oo , 



Indeed, assuming that these 2 integrals are finite, Ito's formula shows that we have for 
some real Brownian motion W" : 



T 2 = T 2 + 2cr W" [ J\t? + 1 - J) dt 



+ 4a*[ S T?dt + 2[ S (r t -f)T t b 2 t d ± 

J J T'f 



JO 



n ■ 



o n 



where 7. , 7' (since |r*t — = 0(6 



-l/2x 



o(a t 



1/2, 



by (iv) and (ii)) are bounded converging 



processes on A . Whence lim T s = 00 almost surely on A , which with (Hi) above implies 
that A must be negligible. 

(vi) End of the proof of the convergence of r s to 00 on {D = 00} . 

By Schwarz inequality, the first bound in (y) above will follow from / T?b t dt < 00 



and from / (r t — ^) 2 b^ dt < 00 . Now these two terms appear in the It 6 expression for 
J 



Zl:=(r s -f)T s b s : 



Z\ = Zl + M\ + 4 



■ + 



dt 



3R\ 2 h 3 Zl — 3. 
2 ! U t „4 2 
' t 



(r t -f)b t , ■ 



dt 

n 



Z]dt+ / T 2 b t dt + / (r t - 
Jo Jo 

with a local martingale M. 1 having quadratic variation : 

(M% = a 2 f\r t - f f b\ x (l - f t + [4 + r t V] T 2 ) eft . 
J 

Note that by Schwarz inequality, (Hi) above implies that / \Z} \ dt < 00 , and then 

Jo 



that 



/* 

Jo 



8 + t 



Z^cft is bounded and converges, almost surely on A , as s — > 00 . 



Using the first assertion of (if), observe that lim 
if (M 1 )^ = 00 . Note moreover that in this case 

dt 

r 2 



[ S (r t -f)ni$+[ S T t \dt 
jo * J 



eft 



00 
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rs (j£ rs 

is also negligible with respect to / (r t - ^) 2 & 3 —7 + / T 2 b t dt . 

Jo rf Jo 

Therefore we deduce from Lemma 3 and from the integrability of t 1— > \Z]\, that 
necessarily (M 1 )^ < 00 , and then that M] has to converge, almost surely on A . 

Hence Z l must have a limit almost surely on A , which must be 0, owing to the 

roo d£ roc 

integrability of Z 1 . This forces clearly / (r t -^) 2 b^ — + T?b t dt to be finite, 

J T J 

almost surely on A , showing the first bound in (v) above. 

Finally, the integrability of T 2 b t and the second convergence of (iv) imply the second 

rOO 

bound in (v) above : / T 4 dt < 00 almost surely on A . 

Jo 

This concludes the proof of the first assertion in Theorem 2. 

2) r s — > R and r s — > 00 occur both with positive probability, from any initial condition. 

Let us use the support theorem of Stroock and Varadhan (see for example ([I-W], 
Theorem VI. 8.1)) to show that the diffusion (r.,b.,T.) of Corollary 1 is irreducible. Since 
we can decompose further the equations given in Proposition 1 for (r.,b.,T.) , using a 
standard Brownian motion (it;., (3., 7.) G M 3 , as follows : 

2 cr 2 r 2 

dr s = T s ds , db s = ab s dw s + a r s df3 s + ^-b s ds + — -r^ (is , 

2 2 6 s 

, U2 p 

rfT s = a T s + a y/l-% d ls + ^T s ds + (r s -\R)^ds-—ds , 

we see that trajectories moving the coordinate b. without changing the others, and tra- 
jectories moving the coordinate T. without changing the others, belong to the support 
of (r.,b.,T.) . Moreover we see from Section 4 that there are timelike geodesies, and then 
trajectories in the support, which link r to r' , and then considering the velocities also, 
which link say (r,b",T") to (r' ,b" ,T m ) . So, for given (r,b,T) and (r',b',T') in the 
state space, we can, within the support of (r.,b.,T.) , move (r, b, T) to (r,b",T"), then 
(r,b",T") to if, b", T'") , and finally move (r',b",T"') to (r', 6', T") , thereby showing the 
irreducibility of (r.,b.,T.) . 

This implies that it is enough to show that for large enough r , T , the convergence to 
00 occurs with probability > 1/2 , and that for ro close enough from R and To negative 
enough, the convergence to R occurs with probability > 1/2 as well. Now this can be done 
by a classical supermartingale argument using the process 1/|T S | , stopped at some hitting 
time. Indeed we see from Proposition 1 that 

1 f s /„2 m o /- on RTt . /„ „ „s bfTt \ dt 



+ f (£z?-^(l-S)-P+(2r,-3J2) 



|T S | Jo V 2 * v r * 7 2r 2 v n 4 ^ |T*| 3 

is a local martingale. 
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Take first r > 3R/2 , T > 4 + jfa , and r := inf [s > I T s = 2 + } : r s increases 
on {0 < s < t} and then we see that l/|T sAr | is a supermartingale, which implies that 

( 2 + ihr'nr < oo) < liminf E(^l {r<oo} ) < liminf E(^) < E(^) < (4 + ^ , 
and then that P( lim r s = +oo) > P(r = oo) > 1/2 . 

Conversely take r < 3R/2, T < —2, and r' := inf js > T s = — v^j: r s 
decreases on {0 < s < r'} and then we see that l/|T sAr /| is a supermartingale, which 
implies that 

2 -i/2 p / r / < oo) < i imin f E(p^ l{r'<oo|) < hminf E( W r^- r ) < E(^) < 1/2, 

and then that F(D < oo) > P(r' = oo) > 1/v^- 

This concludes the proof of the second assertion in Theorem 2. 

3) Existence of an asymptotic direction for the relativistic diffusion, on {D = oo}. 

We want to generalize the observation made in Section 2.1 for R = 0, see Remark 1. 
Recall from Lemma 2 that it does not matter for this asymptotic behaviour whether we 
consider the trajectories as function of s or of t(s) (id est as viewed from a fixed point). 

We shall use Remark 1 and proceed by comparison between the flat Minkowski case 
R = and the Schwarzschild case R > 0. Let us split this proof into four parts. 



f°° a t f°° Ut 

(i) We have / dt < oo and / — dt < oo , almost surely on {D = oo}. 
Jo rf Jo r t 



r°° Ut 

r -t Jo n 

We know from 1) above that r s — > oo almost surely on {D = oo}. 



The very beginning of this proof remains valid : Using (l,i) again, we have almost surely 

roo ^2 «oo ^2 «oo ^2 

/ — ^ dt and / (1 — f^) 2 — ^ dt finite, whence / — ^ dt finite, and then, since 
Jo a t r 2 Jo V 2rt > a t r? Jo a t rf 

a T 2 b 2 1 , r°° a t , n 

— < — - H H - , also / -^dt finite, almost surely on \p = oo}. 

r 2 ar 2 ar A ar 2 Jo rf 

U b a 

Now by the unit pseudo-norm relation, we have — = — < , whence 

r H r 2 i _ R 



r°° U t 

/ — dt finite, almost surely on {D = oo}. 

JO T t 



Ut 
r t 

(ii) The perturbation of the Christoffel symbols due to R is 0(r~ 2 ) . 

Recall from the beginning of Section 2.3 the values of the Christoffel symbols T l jk . 
Denote by Y^ k the difference between these symbols and their analogues for R = 0, which is 
a tensor, has only five non- vanishing components in spherical coordinates, and then is easily 
computed in Euclidian coordinates (x\ = r sin tp cos ip ; X2 = r sin tp sin ip ; x% = r cos ip) : 
we find 

~ . = da^ / dr Or dip dip x 

xj,x k dr \dxidx k rr dxjdxk w dxjdxk 
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since 



dip 



r v 2r(r — R) r r oxj ox k 0Xj0x k ' 

:)r < 1/r and < 1/ (r sin ip) . The same is valid directly for the remaining 

components T^ t and f 1 * t . 

(m) 77ie parallel transport converges, almost surely on {D = 00} . 

Recall from Theorem 1 (in Section 2.2) that the inverse parallel transport £ (s) along 
the C l curve (£ s / | < s' < s) satisfies 

so that, using (n) above and \T S \ < a s : 

t(s) = £o(r; 2 x\i v \)dv = j°0(\i v \ + \r v \+r v \6 v \) r fdv = J°0(2^ + ^)dv. 

Hence, using (i) above, we conclude that the parallel transport (as its inverse £ (s)) admits 
a finite limit as s — > 00 , almost surely on {D = 00}. 

(if) End 0/ i/ie proof. 

Recall from Theorem 1 (in Section 2.2) that the continuous process ( Cs = £ ( s ) is ) is 
a hyperbolic Brownian motion on the hyperbolic space T^M., isometric to H 3 . 

Now, according to Section 2.1, where merely a s = p° , we know that Cs/°s converges 
almost surely as s — > 00 towards (1,1,^) (in coordinates (t,r,0)), for some random 
9oo € S 2 . Using (m) above, we deduce that £, s /a s converges almost surely as s — > 00 
towards (1, 1, ^00), fo r some random ^ G S 2 . This means also that the velocity dZ t /dt of 
the trajectory Z. := (t s \— > (r s ,^ s )) converges almost surely towards (1,^), 1 being here 
the velocity of light. Merely integrating this, as in the proof of Remark 1, we get finally 
the generalization of Remark 1 to the relativistic diffusion. This ends the whole proof of 
Theorem 2. o 

Note moreover that the part (3,(iv)) of the above proof shows that, since the hyperbolic 
Brownian motion does not explode, there is no explosion at the level of the fibre. This 
is a general fact for relativistic diffusions. Finally if the radius r s could explode within 
some finite proper time, T s , and by the unit pseudo-norm relation the energy a s , would 
explode as well ; but this is clearly forbidden by the simple stochastic differential equation 
governing a s . This proves the non-explosion of the Schwarzschild relativistic diffusion, 
quoted directly after the statement of Theorem 2. 

3.2 Proof of Theorem 3 

(i) The first sentence is clear from Theorem 2 when the diffusion starts above Sq , 
except the finiteness of D', proved in (ii) below. This is the same when the diffusion 
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starts above — So , Theorem 2 being valid as well in this very similar case (it is sufficient 
to change the signs of a and t). The other cases are reviewed in Remark 6. So we just 
have to establish the assertions (1), (2), (3) of the statement, assuming D finite. 

(ii) Let us prove (1) first. Since T D < and T s cannot vanish in the region {r < i?}, 
r s must decrease and then converge to some g G [0, R[ , with limsupT. < —\J^ — 1 < , 
which in turn forces g — , hit within a finite proper time D' , and lim s/ *£>/ T s = — oo . 

Only the upper bound for D' remains to be proved. It is a consequence of the pseudo- 
norm equation, which implies T 2 > y — 1 and then T s < —^J ^ — 1 on {D < s < D'} . 
Indeed, consider 



g(r) := f Arcsin ^r(R-r) -^r(R-r) if 0<r<|, 



and g(r) := f (tt - Arcsin [^y/r{R - r) ]) - ^Jr(R-r) if f < r < R . 

We have 5 '(r) = (f - l)" 1 / 2 , <?"(r) = ^ x (f - I)- 3 / 2 , ^(0) = </(0) = , 

(?(/?) = I i? , = +°° • Hence g'(r s ) x T s < -1 implies by integration on [D, s] , 

for D < s < D' : ^(r s ) < f R + D - s ; taking s = D' , this proves (1). 

(iii) Let us then prove the non-explosion of log b. at D'. For that, let us apply the 
comparison theorem (see ([I-W], VI, th 4.1)) to 6[d,d'[ : we get so a real diffusion process 
P. solving 

p s = bD + a f D y/ff + R? dw s + "4 f D {p s +R 2 )j g , 
and such that almost surely sup b s < sup f3 s . Indeed the ratio of the drift coeffi- 

D<s<D' ' D<s<D' 

cient and of the squared diffusion coefficient of b. is maximal for r s = . Now there is a 
real Brownian motion B. such that 

fe+ ,,^-^[ M{s |/; _^_ >4 , }] , 

so that j3. cannot diverge at a finite time ; indeed it is immediately seen that 

I —j = 21og(S /.B s ) + 2 / — — must almost surely diverge as s approaches the hitting 
J lif Jo B t 

time of by B . This proves that b. almost surely cannot explode at D' , and thus must 
be continuous on [0, D'\ . 

Moreover, since the differential equation governing b. can be written 

r r s / T 2 \ 1 

\og(b s /b ) = aW / + +a 2 s, for some standard Brownian motion W. , we see 

l jo v bi> J 

rD' r 2 

that bo> = would imply / -y dt — 00 , and then limsup 6 S = +00 , a contradiction. 

Jo s /D> 

(iv) Let us now prove the statement (3). For that, let us write again the unit pseudo- 
norm relation : near D' it writes 
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T s r 3 J 2 = -^{R - r s )b 2 + Rr 2 + (a 2 s - l)r 3 —+ -b D ,y/R . So that for any e > we 

have: —bu'\fR{^ + s) < T s r 3 J 2 < — 6 D /V r R(l — e) , for s sufficiently close to £>'. 
Integrating this, we get immediately, for s sufficiently close to £>' : 

-bD>y/R(l + e)(D'-s) < -frf 2 < -b D ,^R (1 - e)(D' - s) , 

which means the first equivalent in the statement (3). The second equivalent follows at 
once by using again the unit pseudo-norm relation. This proves (3). 

As a consequence, we deduce at once that \9 S \ = U s = b s /r 2 s is integrable near D 1 ', 
which implies the convergence of 9 S in S 2 as s / D' '. 

(v) End of the proof of the statement (2). 

To establish the non-explosion of a s , let us rewrite its stochastic differential equation 
with two independent real standard Brownian motions w,/3 : 

da s = aa s dw s + a(f-- 1) df3 s + ^a s ds, and consider X s := a s x e -<"».-D-o a (s-D) _ 
We have dX s = e -°™s-D-° 2 (s-D) a - 1) d(3 s , whence for D < s < D' : 

tts = e ^s- D +«Hs-D) ^ +a J\-a Wt _ D -„^t-D) _ x) ^ 

for some real standard Brownian motion VK. The equivalent seen above for r s near D' 
shows the almost sure convergence of the integral fp , and then by the above formula, 
of a s as S D'. 

The same equivalent again shows the almost sure integrability of ii~ near D' : indeed 
it; = (a s + T s )/{1 - R/r s ) ~ [2 5 4 D ,/(5 R 3 (D' - s))]" 1 / 5 , 
which finally proves the convergence of u~ as /* D' . Likewise for uf . 

— * — * 

Finally it remains to show the convergence of b s , or equivalently of b s /b s = 9 S A0 s /U s , 
since we already saw above the convergence of b s in M.* + . Since we also saw the convergence 
of 9 S in § 2 , it remains to get the convergence of 9 S /U S . Now this is the last assertion of 
Proposition 2. o 

3.3 Proof of Proposition 2 



We have 



d(?f) = 2 r s b- 1 T s 9 S + r 2 b; 1 d9 s - (r 2 bf db s ) 9 S + r 2 s bf (db s ) 9 S - r 2 bf (db s , d9 s 



To perform the computations, let us use the basis (u,v,k) of R 3 defined by: 
u = (cos if), sin ip,0) ; v — (— sin ip, cos tp, 0) ; = (0, 0, 1) , 
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so that 9 = u sin tp + k cos tp , and 

9 — (u cos <p — k sin ip) tp + v ip sin ip ; 9 A 9 = v tp — (u cos ip — k sin ip) ip sin tp . 
Then d# = (u cos </? — sin (p) dtp + v sin tp dip — tp 2 9 + 2 v tpip cos ip — u ip 2 sin ip , 
whence (d#, d&) = (u cos y? — k sin </?) (d</?, d&) + v sin 99 (d^, db) 

= ((u cos (/? — k sin + v sin </? ^ <x 2 b^ 1 (b 2 + r 2 ) = b^ 1 (db) 9 , 
so that the last two terms in the expression of d(j%) above cancel. Hence 

d(^) = 2 £ T s S - r 2 bp 2 (a dM b s + b s ds + ^ ds) £ s 

+ C/ s _1 (« s cos tp s — k sin </? s ) (er dMf + ^ tp s ds — ^- T s tp s ds + sin </? s cos y2 s ?/> 2 ds) 
+ C/ s _1 v s sin y? s (o- dM^ + ^fip s ds - ^T s ip s ds - 2 cotg ^ s rf s) 
+ Upp 1 {2 v s ip s ips cos ip s - ip 2 s 6 s - u s tp 2 sin tpfj ds 



( — 9 S dM h s + (u s cos ip s — k sin </? s ) 6 S dMf + t> s 6 S sin 



2 4 



6? 

1 7 .j # s ds — L7 S _1 V^s ^« ds + C/ s _1 sin </? s ((u s cosy? s — A; sin tp s ) cos ip s — u,fjip 2 ds . 
2 b s 

Observe now that the definition of b = r 2 U implies 

dM b s = r 2 Uf 1 (tp s dMf + sin 2 tp s 4 dMf) . 
Therefore, expressing all in the basis (9 S , 9 S /U S , 9 S A 9 S /U S ) , we get 

d(^) = a Uf 2 sin tp s (y 8 dMf - ^ dMf) 9 s A^-^-9 s ds- U s 9 S ds , 

that is to say the formula of the statement, with d(3 s := r s C/" 1 sin (p s (tp s dMf — ip s dMf) . 

Now it is straightforward to verify (from the covariation matrix given about the begin- 
ning of Section 2.4.1, before Theorem 3) that 

(d(3 s ) = ds , (d(3 s , dM h s ) = , (d(3 s , dMj) = , 

which shows that indeed (5. is a standard Brownian motion and is independent from 
(r-,b.,T.). 

To establish the last assertion of the statement, we deduce from this expression for 
d^ff) that the bounded variation part of 9 S /U S is not larger than b s rj 2 + cr 2 b~ 2 r 2 , which 
(as r~ 2 , recall Theorem (2,3)) is almost surely integrable on [0, D'\ , while its martingale 
part has quadratic variation not larger than a 2 bp 2 r 2 , and thus is almost surely integrable 
on [0, D'] as well by Theorem (2,2). o 
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3.4 Proof of Proposition 3 



1) To establish this, we need a system of stochastic differential equations relative 
to the whole relativistic diffusion (£ s ,£ s ), and not only to its projection (r s ,b s ,T s ) on 
the coordinates (r, b, T). Recall that the whole relativistic diffusion (£ s ,£ s ) lives in a 
7-dimensional space. For our purpose the following system of coordinates is convenient : 
(r, 9, u~, a, T, 9) E R + x§ 2 xlxlxlx T § 2 . Recall that 

b = r 2 \9\ = r 2 U = r 2 \j<Jp 2 + ?p 2 sin 2 <p . 

Now we have the following simplification : we can recover the coordinate u~ from the 
other ones by integration, since its value is fixed at the origin. Hence, to recover the whole 
relativistic diffusion (£ s ,£ s ), it is sufficient to get the reduced diffusion (r s ,a s ,b s ,T s ), and 
to recover (9 S ,9 S ) conditionally from (r s ,a s ,b s ,T s ). 

2) Let us deal with the reduced diffusion (r s , a s , b s ) : we have to show first that for any 
(cto,&o) ^ Kxl^. there exists a unique (in law) diffusion process (r s ,a s ,b s ) G [0,fl]xRxR^ 
defined up to proper time D" : = inf{s | r s — R} , starting with initial condition (0, a Q , b ) , 
and having infinitesimal generator (see Proposition f ) 

where T = T(r,a,b) := \ja 2 + (f — 1 )(^ + 1) is chosen positive. 

Clearly r s must increase strictly as s / D" , so that s can be expressed as random 
continuous strictly increasing function of r e [0, R] : s = s(r) . Let us set a r := a s M and 
b r := b s M . In other words, we consider here the radial coordinate r as an alternative time 
coordinate. Then (a r ,b r ) has to be a time inhomogeneous diffusion on R x R+ started 
from (ao,&o) and with infinitesimal generator 

£ := ( (S 2 - f + ^ + (b 2 + r 2 ) + Tab -^L + 35 + (36 + ■£) . 

Conversely, given such a diffusion, the inverse time change yields a diffusion with gen- 
erator Q" . Therefore it is enough to prove existence and uniqueness for the (/-diffusion 
stopped at r = R. Now it is well known that a sufficient condition is locally boundedness 
and continuity of the coefficients of the associated stochastic differential equation, together 
with a local Lipschitz condition on these coefficients with respect to (a, 6) . 

Now the term r 2 /b causes no trouble since bo ^ and since the proof of Theorem 
(2,2) insures that b r stays in R*j_ , and then we observe that (a 2 — I + ^) jT goes to as 
r \ , and that T" 1 / 2 , T -1 , and their derivatives with respect to (a, b) , stay bounded 
as well ; indeed (for example) T~ 3//2 | g | = 2 T^' 2 (f - 1) ^ < 2 T' 1 ' 2 /b is bounded. 

3) It remains to prove that we can recover (9 S , 9 S ) conditionally from (r s , b s , T s ), once 
(9q, #o) is fixed, in a unique way. Now Proposition 2 displays the equation we have to 
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solve. To solve it, let us complete the equation of Proposition 2 into a linear system in the 
variables V s := (9 S , (6 S /U S ), 9 S A (9 S /U s )^j, by adjoining the equation d6 s = U s (6 S /U S ) ds , 
and the following one (immediately deduced from Proposition 2) : 



17, 



So we get a linear differential system : dV s = V s dA s , where the matrix- valued differential 
dA s is given by : 



dA s :-- 



~(b s /r 2 s ) ds 

(b s /r 2 )ds -{a 2 r 2 s /2b 2 s )ds -(ar s /b s ) df3 s 
(ar s /b s )dp s -{a 2 r 2 s /2b 2 )ds, 



Note that this differential system can be equivalently written in the Stratonovitch form 
dV s — V s o dA s , where the matrix- valued differential dA s is given by : 



dA s :-- 



-(b s /r 2 )ds 

(b s /r 2 )ds -(ar s /b s )dp s 

((rr s /b s ) d[3 s 



so that any solution takes its values in the rotation group. 

Let us solve this linear equation by means of the following series : 



V s = V D i ( 1 + ^fc(s)) , where for each k Jk(s) := / 

keW {D ' 



<si<..<s k <s} 



To justify that, let us choose on the space of (3, 3)-matrices the Euclidian operator 
norm, and fix C — C(u) > 1, measurable with respect to (r.,b.,T.), such that for 
D' < s < min{D", D' + 1} we have : 



b s /r 2 <C(s- DT 4/5 , °r s /b s < C (s - D'f 5 , a 2 r 2 /b 2 s < C ; 



< 



(5C) 2 



A;! 



x (s-D') 2k /\ 



this is possible by Theorem 2. Let us suppose that E \\J k (s) u2 
where T° denotes the a- field generated by the reduced diffusion (r., b.,T.) . Then we have 

2 



E 



JD' 



E 



Jk+i(s 



= E 



dA s , 



i\ 1 1 2 



f Jk(s') 

JD' 

'fLf ds ' + 2 f k^ds'x [ S E\\\J k (s') 

JD' TZ, JD' L 



b,' 



ds' 



+2(s-D') x [ S E\\\J k (s') 

JD' L 



2J2, 



,o r 



s \2 J J 



2b 2 



ds' 



<2C 2 



(5C) ' 
k\ 



2k r pg 



f (</ _ jD ')(2fc+4)/5 ds , + 5(s _ f S (s / _ £/)(2fc-4)/5 ^ 

JD' JD' 

+ {s-D') f S (s'-D') 2k / 5 ds' 

JD' 
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— 2 C 



k\ 



2k + 9 



+ 5 



2k + 1 



+ 



2k + 5 



(r.n\2{k+l) 



converges, this shows by induction that the series Jfc(s) 



Since the series > 

* (*0 1/2 

almost surely converges in L 2 , conditionally with respect to the reduced diffusion (r., b.,T.) . 

As to the unicity of this solution, note that any other solution V' s must satisfy for any 
nGf: 



V. = V D , (l + ± J k (s)) + / 

k=l J ^ u 



' {£>' <s <si<..<s n <s} 



V' x dA S0 x dA Sl x .. x dA Sn , 



so that, using that ||V^ o || = ||Vd'|| i s constant, we get (with the L 2 -norms understood 
conditionally as above) : 



\v;-v s \\ 2 < 



v; o dA so ..dA Sn +\\ y: j k (s) < nj n+1 ( s )ii 2 +x: ii j k ( S )h 



{D'<s < S1 <..<s n <s} " u ' '• 2 II k>n 2 ^ 



(5C) k 

< 2^2 m/2 for any n E N* , showing that V' s = V s almost surely. 

k>n ( k -) 

This concludes the whole proof of Proposition 3. o 



3.5 Proof of the second assertion of Theorem 4: = oo a. s. 

For any e e]0, 1[, set A £ := { sup b < e^ 1 } ■ By continuity of b, we have 

P({Doo < oo} fl A°) < e' , for any fixed e' > and small enough fixed e . 

For any nGN, set r' n := inf{s > D 3n+2 | T s = 0} and r n := inf{s > D 3n+3 | r s = R/2 } . 

Note that D 3n+2 < r' n < D 3n+3 < r n < D 3n+4 on {D^ < oo}, and £ (r„ - < . 

n 

Moreover r n — r' n = r o @ T ™ , where r denotes the hitting time of {r = R/2}. 

Setting A' £ := {r < ce} fl { sup b < e^ 1 \ (for some constant c = c(a, R) > ), we have 



P(D 00 < oo) - e' < P({,D 00 < oo} n A £ ) < P(^(U S x r) o 6 T » < oo) 

n 

P(lim n inf(e<)- 1 (^)) <£ hmP( f] (0^)^(4)) 



< 
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by the strong Markov property. Hence we see that this proof is achieved if we show that 
lim E(P e (A'J x n (Iaj o @ Tra ) = 0, for any n. Now this follows immediately, by 

p—*oo \ T p m=n £ ' 

induction, from the following lemma, o 

Lemma 4 For small enough e > , and for any initial condition r > R and b > 0, 
provided T = ; A' £ being as in the proof of Theorem 4 above, we have P(^) < 1/2 . 

Proof Let us write again the stochastic equation of T s , under the following form, for 
some real standard Brownian motion W : 



T s = aW 



Consider r := inf{s | r s = R/2} and o' e := inf {s | |T S | = e *} , and some constant q . 
On the event A" := \ max |W|([0, (e- 2 + l)g]) < (2ae)- 1 \ n { sup 6 < s" 1 } , we have 

for < s < minis, a' r} : r s = r + T t dt > R — e~ 1 s , and by the equation of T s : 

Jo 

-h ~ "4 ? - 3 (2/i?) 3 ^ 2 q ~ (2/R) q<T s <j- £ + 3 4q+ (2/i<>) V 2 q , 
1 1 

whence \T S \ < — + — (3<7 2 + 4 iT*£ + 48 RT 3 e -1 ) g . 

Hence g < min{o^, r} on A" £ if g < (3<r 2 + 4i?~ 1 e + 48 i?~ 3 e _1 ) 1 and q<Re/2. 



Thus there exists a constant c = c(<7, i?) > such that min{o"g, r} > q = ce on A" . 
Therefore { max |W|([0, (e~ 2 + l)g]) < (2(T£)- 1 } n A' e = , and then 

P(A^) < P[max|W|([0, (e' 2 + l)g]) > ^ae)" 1 ] = P[max |W|([0, g]) > (2oVe 2 + 1 ) _1 

< 2P[maxH/([0,g]) > (2a Ve 2 + 1 = 2¥[\W 1 \ > (2a Ve 2 + 1 / y/ce' 

= 4f e -*v* *l < ^L^e-y^. o 

Jl/(2a^/ce(l+e 2 )) \/2lI V2n 

3.6 Proof of Theorem 5 

The proof of this theorem will be split into 8 parts and a series of lemmas. 

1) Let us begin by the dichotomy of the first assertion: there is no other possibility 
than the obvious one : r. goes to infinity, and the confinement exhibited here : r. remains 
endlessly bounded. 

Lemma 5 Almost surely, if r. does not go to infinity, then it is bounded. 
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Proof Consider the double sequence of hitting times {A n , AJ, | n 6 N} defined by A = 
Aq = , A' n : = inf{s > A n \r s < R} , A n := inf{s > A' ri _ 1 | r s > e n and T s = 0}, and 
set E n := {A n < oo}, for n G N*. 

By Theorem (2,1) the event fl n E n contains all trajectories such that r s does not go 
to infinity, but is unbounded. Thus we want precisely to prove that P(n n E n ) = . 

Let us apply the comparison theorem (see for example ([I-W], Theorem 4.1)) : there 
exist comparison processes and X~ on the same probability space, such that almost 
surely on E n : 

max T s > max Xt , Xt = and min XL > min X~ , X7 = , 

A„< S <A„+1 - A„<s<A n +l A " A n <s<A„+l ~ A„<s<A n +l 

where X^ and X~ are real diffusions given by their stochastic differential equations, 
which are deduced from the equation governing T s by bringing down the ratio between 
the drift term and the diffusion coefficient, and by bringing the diffusion coefficient down 
(in the X + case) or up (in the X~ case) . Recall that 

dT s = a^T* + l-f- dw s + ^T s ds + (r s -f)^ds-^ds. 

Thus, as long as r s > 3R/2 , we can take | < T 2 + | < T 2 + 1 - f < T 2 + 1 and then : 

dX+ = ^dw s + 3|i (X+ A0)ds-^ds, 

and 



dx; = a y/(X 7 y + ldw a + (x; A 0) - ds . 

Now, for A := 4 + , fix e := Pf max Xt > A) , which is > , > such 

that Pf min X~ < -N) < e/2, and n such that r An > N + 3R. Note that the 

event E' := j max Xf > A} fl j min Xr > —N\ has probability > e/2. 
n U„<s<A„+l s J La„< s <a„+i s - ) 

Setting A" := inf{s > A n \r s < 3R/2} and applying the comparison theorem, we get : 
min T s > min X~ > -N on E n C\E' 

A„<s<(A n +l)AA^ A„<s<A n +l 

and then 

min r s > r A - N > 3R on E n C\E' 

A„< S <(A„+1)AA» - " 

showing that A^ > 1 on _E n fl E 1 ^ . Hence, applying the comparison theorem again, we 
get also : 

max T s > max Xj~ > A on E n fl E' . 

A n <s<A„+l ~ A n <s<A„+l 

Thus, using the strong Markov property, we find that for large enough n : 
P((3 s > A n ) T s >4+j^ and r s > 3R/2 / E n ) > e/2 . 
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Let us use now the proof of Theorem (2,2), where we proved that for r > 3R/2 and 
Tq>A + , then P( lim r s = oo) > 1/2 , together with the Markov property : we obtain 

f(E n+1 1 E^j < p(A'„ < oo J E^j < 1 - e/4 , for n larger than some n . 
Therefore we get F(E no+k+1 ) < (1 - e/A) k and then F(n n E n ) = , as wanted, o 

2) The first case in the theorem was already handled in Lemma 2 and Theorem 2. The 
irreducibility of the relativistic diffusion is clear from Theorem (2,2). 

Let us now focus on the second case in the theorem, supposing therefore that r. is 
bounded. Set r M ■= inf{s | r s = M}, for M of the form M = k x 3R/2 , with k G N*. 

The proof is divided in several distinct lemmas. We shall always let the relativistic 
diffusion start at proper time D_i = from level {r = R} with Tq > . We let the 
hitting times Dj be as in Section 2.4.3. Fix also some e in ]0, |[. 

3) Estimates related to b. 

Let us begin by proving that (when r. is bounded) b. has to go to infinity, which by 
means of the Markov property will allow then to consider only large enough b . 

Lemma 6 Almost surely, if r s is bounded, then lim b s e~^ l ~ £ ^ a2s = +oo . 
Proof Note that {sup r s < oo} = (J {tm = oo}. 

s>0 M 

Let us recall the logarithmic form of the stochastic differential equation of b : 

\og(b s e _ ( 1 ~ £ ) cr2s ) = log&o + £o~ 2 s + aW s + / ^| dt (for some Brownian motion W). 

Jo b t 

A straightforward consequence is that limsup(6 s e~( 1_e )°" 2s ) = +oo almost surely, so that 

s— >oo 

the stopping time t n := inf {s > | b s e _ ( 1_£ ) CT s > n} is finite for any n G N . 

Then let us write the equation governing b. in the following linearized form (for some 
other Brownian motion W) : 

hse -{i-e)^s = eaWs+ e^ s ^ + a w[j\- 2 ° w ^ 2 ° 2t r 2 t dt] + j* e- awt - a2t ^dt). 

It is clear from this expression that for any fixed M the probability of the asymptotic 
event A M := {t m = oo} fl {b s e^^ £ ^ (J s does not go to oo} can be made arbitrarily close 
to by choosing b large enough. Using that A M = Q^(A M ), we apply the strong 
Markov property at the sequence of stopping times {t n } to conclude that F(A M ) = , 
which yields the result, o 

We estimate then the increase of b , when started from some large value b . 
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Lemma 7 Fix e > 0, and recall that tm '■= inf {s \ r s = M}. Then there exists a lower 
bound b(M, e) such that for b > b(M, e) we have : 

F(bl {£/2) e^ 2s <b s < bl +(£/2) e^ 2s for all s < r M ) > 1 - 2 b'/ 12 >l-e. 

Proof For any q g]0, 1[ , set v q := inf{s | b s e -( l ~ £ ) a2s < qb ] g]0, oo] . 
Recall that the equation governing b. can be expressed in the following form : 
r f s / i 

\og(b s /b ) = aW[^ (l + -±)dt\+o- 2 s , for some standard Brownian motion W . 
Hence for < s < u„ A r M we have : 



whence 



log(6 a e-^-^V&o) > ea 2 s + a min W[s, (1 + (M/qb f)s] , 



PK < t m ) < P((3 s) e s + min W[s, (1 + (M/qb f)s] < log q) 



Taking q = b^ e/2 , this yields : ¥(u q < t m ) < ~ £2 ^ 1+mH ^ < ^ 2 / 2 f or ^ > b(R, e ) . 

This proves the lower control of the statement. The upper control is obtained exactly in 
the same way. o 

4) We prove now that q := limsup r s cannot be strictly between 3R/2 and +oo . 

Lemma 8 Almost surely, if r s is bounded, then g := limsup r s < 3R/2 . 

Proof Let us proceed somewhat as for Lemma 5. Fix e e]0, 1[, M > 3e + 3R, and 
consider the double sequence of hitting times {A n , | n G N} defined by A = A' = , 

A n := inf{s > A' n _ 1 1 r s > 2e + 3R/2 and T s = and b s > e 2n }, 

A' n := inf{s > A n \ r s < e + 3R/2 or b s < e n }. 

Recall that t m ■= inf{s > | r s > M} , and set E n := {A n < r M }, for n G N*. 

By Theorem (2,1) and by Lemma 6, the event l~l n E n contains all trajectories such that 
r. is bounded by M and limsup r s > 2e + 3R/2 . Thus we want precisely to prove that 

s^oo 

P(n n E n ) = o. 

Let us apply the comparison theorem (see for example ([I-W], Theorem 4.1)) : there 
exist comparison processes , X~ and Y~ on the same probability space, such that 
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almost surely on E n : min b s > min Y Q , Y" A = 6a„ > e 2n , and 

max TL > max , Xa~ = and min T s > min XZ , X7 = , 

A n <s<A n +l A n <s<A n +l n A n <s<A n +l 0<s<l 

where and X~ are real diffusions given by their stochastic differential equations, 
which are deduced from the equation governing T s by bringing down the ratio between 
the drift term and the diffusion coefficient, and by bringing the diffusion coefficient down 
(in the X + case) or up (in the X~ case), and similarly for b s , Y~ . Recall that 

2 2 

db s = a \Jb 2 + r 2 s dw' s + ^-b s ds + — — - ds , and 

Z Of, 



bi , R 



dT s = a JT? + 1 - - dw s + ^T s ds + (r s - ^) ^- ds -ds. 



2 2^"" 2 r 2 



Thus, on [A n , A' n ] we can use b 2 + r 2 < b 2 + M 2 and T 2 + § < T 2 + 1 - f < T 2 + 1 , 
and then, for n such that eM" 4 e 2n — ^ > 3 e n : 



dy s " = a ^{Y-) 2 + M 2 ^ + 3 4 Y- ds , 
dX+ = a y/(X+) 2 + \ dw s + ^-{Xj AO) ds + e n ds, 

and 

dX; = <r\l{X-) 2 + 1 du> s + 5f (Xf A 0) ds + 3 e n ds . 
Now Pf min Y~ > e n ) , Pf min X" > -s) , and (for A := 4 + -A ) 

U„<s<A„+l s / U„<s<A„+l s / Ha 

P( max Xj~ > A) are (for large enough n) arbitrary near from 1, so that the event 

<s<A n +l / 

F n := i min Y~ > e n \ f] \ min Xr > — e \ fl i max Xt > A\ has proba- 

L A n <s<A n +l J L A„<s<A„+l s J L A„<s<A„+l s J 

bility arbitrary near from 1. 

Hence, applying the comparison theorem, we get : 

min b s > min Y~ > e n on F n fl E n , 

A„<s< (A„+1)aA^At A// A„<s<A„+1 

min T s > min X~ > — e on F n fl E n , 

A„<s<(A„+1)aA^Atm A n <s<A n +l 



and then 



min r An+s > r An - e > e + 3R/2 on F„n£ n , 

A„<s< (A„+l)AAJ l Ar M 



showing that A' n > A n + 1 on F n fl {tm > A' n } fl £■„ . Hence, applying the comparison 
theorem again, we get also : 

max T s > max X+ > A on F n n |r M >A!,}n£„. 

A„<s<A„+l - A„< S <A„+1 s L - nJ 
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Thus, using the strong Markov property, we find that for large enough n : 
P(r M < A; or [(3 s G [A n , A^[) T s > 4 + ^ and r s > 3R/2] / E n ) > 1 - e. 

Let us use now the proof of Theorem (2,2), where we proved that for r > 3R/2 
and T > 4 + , then P(r s increases towards oo) > 1/2, together with the Markov 
property : we obtain 

p(tm < A' n I E^j > 1/2 — e , for n larger than some n . 

Therefore, noting that E n+1 C {r M > A' n }, we get F(E no+k+1 ) < [1/2 + e] h and then 
P(fl„ = , as wanted, o 



5) Lower estimates related to D 3n+ 2 

The aim of the following lemma is to get a lower bound (in an optimal way, owing to 
Remark 12 below) on the duration of an excursion in the hole, ouside an event of small 
probability. The idea is that when T. vanishes, with an acceleration of order b 2 it takes a 
lapse of time of order 1/6. to make a non-infinitesimal move. 

Lemma 9 Suppose r = R , T > , and b large enough. 

Then we have ¥(d 1 >^At m )>1- e -V(64<x 2 c) ? y or c ;= R^/t^M). 
v % ' 



Proof Consider D' := inf{s > | T s = 0} < D , D[ : = inf{s > | r s = f } G [D' , D ± ], 
r] : = inf {s > | log(^) = ± log 2}, and A := inf{s > D' \T S = ±^} < D . 



Set E := {D' < ± A r M }, £' := {<7 Vy*[2c/fc ] < §}, and £"' := {a iy*[2c6 /i? 2 ] < f%} 



4RJ ' 



where W and are standard real Brownian motions used to write the stochastic equa- 
tions governing b. and T. respectively, and W* , W* are their maximum processes : 
W*(s) := max |W|([0, s]) for any s > and similarly for IV* . 

^.2 



Recall from the proofs of Lemmas 6 and 7 that \og(b s /b ) = aW 



1 



%)dt 



+ CT-S, 



and (as already used for example in the proof of Lemma 4) that for s > D' : 



T x = aW 



(T/ + 1- Z)dt 



Therefore on E fl E', for < s < ^ A t m A r] we have : 



log(6 s /6o) 



<aW* 



1 + +e<l + e<log2, whence r, > £ A r M 



bo. 



cr 2 c ^ 1 | <j z c 
bo 2 bo 



bo 



and on the other hand, on E fl E' fl E", for -D < s < ^ A r M A A A D' 1 we have : 



ITJ < o-W* 



' _^o_ _i_ -1 _c_ I 3<t 2 bp c I o /\# 64bg c. i 2 _c , jo j 3cr 2 c , 128Mcbp , 2c ^ bp 
>R 2 ^ / 60 J 2 R 60 « 4 bo T 60 4fl T 2fl T iJ 4 Rb R ' 
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since c = i? 3 /(512M) , for b large enough. Since c < R 2 /2 , this implies that 

r s = R + jT T t dt > R - | £ > f , still for < s < £ A r M A A A D[ . 
Hence we find that EdE'dE" C > ^Ar M }, and then also E'dE" C > ^Ar M }. 

P(f > = PfrA^- < f) = P( sup | log(fc) | < log 3, sup |T S | < |, sup |r a -i2| < f). 

Now, as already used in the proof of Lemma 4, we have for b large enough : 
¥[E'} > l-8a ^/c/(nb )e- bo/{WcT2c) > 1-I e -W(i^) 

and P[£"] > 1 - 16a v /c/(7r6 )e- 6o/(64<T2c) > 1 - ± e^A 64 ^ , 
whence finally : 

F(D 1 > A r M ) > P(D; > A t m ) > P(#' H £"') > 1 - e -W(64* 2 c) . o 

Let us apply now the preceding lemma, to get jointly a lower bound on all proper 
hitting times at which an excursion in the hole begins. 

Lemma 10 Fix e > 0, and recall that tm '■= inf{s | r s = M}. Then there exists a lower 
bound b(a, R, M, e) such that for b > b(a, R, M, e) we have : 

P(£>3n-i > r M A log(l + n)) forall neN)>l-e. 

Proof Set 5k '■— e( 1+e )°" 2 D3fe + 2 , for fcsN, and consider the following events, indexed by 
the integer n > : 

B n := {bl^ e^° 2 ° <b s < bl + ^ e ^eWs for aU s < 
A n := {^ 3 fc-i > r M A ((i+i)^ log(l + ( ) fc)) for all integers fc < raj, 

and 

< := {£>3n-i > log(l + (^#) n)} , K := B n fl {r M > L> 3 n-i}- 

Finally let us consider also 

5 := ed-)^- < 6 S < e^ 2 * for all s < r M }. 
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We have 

4*i = {Ds n+ 2 > log(l + (<£$£)(n + 1))} = {5 n > 1 + (G$7?)(n + 1)} 



and 



D {(£>3n + 2 " £>3n-l) <f„-l > "^y} fl {*n-l > 1 + (^f) n}, 

{(£>3n+2 - ^3n-l) #n-l > l+(e/2) } D {^3n+2 ~ ^3n-l > 7 } fl B n , 

On °-D 3n _i 



whence, using that A n+1 = (A' n+1 H A„) U ({-D 3 n+2 > t m } H A n ) : 

A t +1 n ^ D {^ 3 n+2 - > 7-^— OI D 3n+2 > T M } f| A n H ^ . 

Now by Lemma 9 and the strong Markov property we have (denoting by T. the natural 
filtration of the diffusion, and setting a := (64<7 2 c) _1 ) : 



P(L>3n+2 " £>3n-l > 

Hence 



? £>3„-l 



or L> 3 „ +2 > r M ^D3„-i) > E(l - e" 



3n+2 — 



> 



or D 3n+2 > r M }fl^ n5 n) 



> E 

> 



(l_ e -«to>s»-i) X I A 

(l — exp 
= (l — exp 



,.nB' 



> E 



(» 



- 1-6 

e 



l-(e/2) X <T+?K 

- 1 j x l Aj 



«6j- (£/2) ((1 + £ )ca 2 n/6i +(£/2) )^]) x P(A n n B' n ) 



64 



Set £ n := exp -^(a 2 c) ^ (1 + e)^ 6 1+£ , for n > 1 , and e :— e ab ° , so that 

we have clearly, for 6 > b(a,R,e), on one hand ^2 e n < b^ 2 '^ < e, and on the other 

n>0 

hand : 

P(A n+1 n B' n ) > (l - £n ) x p(A n n > P(A» n - ^« , 

or equivalently 

p(A n n n {A n+1 f) <e n . 

Note that 5; D {r M > £>3n-i} (1 B . Therefore 

P(A n+1 n B) = p(A n n {L> 3 n-i > m} hb)+ p(a h+1 n {r M > Asn-i} n £?) 
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= p(A n n{D 3n ^ > T M }nB) +p(A n n{r M > D 3n ^}nB) -p(A n n{r M > v 1 }nBn(A,+i 

> p(A n n {L> 3 n-i > r M } n £?) + p(A n n {r M > d^-i} n b) - p(A n nB' n n (A n+1 ) c 

> P(A n n {D 3 n-i > r M } n b) + p(a„ n {r M > £ 3 n-i} n b) - e n = ¥{A n n B) - e n , 
whence finally by Lemma 7, for 6 > R, £) '■ 

P( fl ^) > P ( 5 ) - E £ n > 1 - 26„ £2/2 - b £(2 - £) >l-e. o 
6) Estimates related to \a.\/b. 

f dt 

We need next to control the integral / — , which occurs in the It 6 expression of the 

J r t bt 

crucial quantity — ^ . The following lemma estimates its contribution due to an excursion 
b s 

in the hole. 

Lemma 11 During any excursion in the hole, id est during any proper time interval 

rDzn+2 dt 71 R 

[D 3n , D 3n+2 ] , we have the following control : — < . 

JD 3n r t mm b 

[D 3n ,D 3n+ 2] 

Proof Firstly, it is sufficient to consider the proper time interval [D 3n , D 3n+1 ] , since the 
estimates are exactly the sames on the other half [D 3n+1 , D 3n+2 ] . Since 

T t n' 2 = -J(R - r t )b 2 t + Rrl + (a? - l)rf < -yfR~=F t b t for D 3n < t < D 3n+1 , we get 



D 3n+1 di ^ rD 3n+1 / y t dr t ^ 1 [ R I d 

D 3n r t ~ Jd 3ti \j R-r t b t ~ min b[D 3n , D 3n+2 ] Jo \jR-r 

ttR 

. o 



2 min b[D 3n ,D 3n+2 ] 



Remark 12 The same argument yields also the following estimate on the duration of 
an excursion in the hole (showing that the estimate from below in Lemma 9 is essentially 
optimal) : 



f Dsn+2 , 2 f R / r Svri? 

D 3n+2 - D 3n = dt < , f x / J — r dr 

JD 3n mm b\D 3n ,D 3n+2 \ Jo 



D 3n min b[D 3n ,D 3n+2 ] Jo y R - r ' 4 min b[D 3n , D 3n+2 ] 

f dt 

We can now deduce the control on the integral / — -~ , which we shall need below. 

J rtbi 



Lemma 12 For any e e]0, \[, there exists a lower bound b(a, R, M,e) such that for 

f TM dt f TM R dt 

bo > b(a, R, M,e) we have : / + / — -k < b^~ 2 , with probability larger than 1—e. 

Jo b t Jo r t b t 
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Proof Recall from Lemma 7 that with large probability we have b s > &o" (£/2) e^-^ 8 

f™ -2 bfr 2 

for s < tm , whence / b t dt < — — — . Moreover, using Lemma 11 we have : 

jo 2(1 — e)a 2 



Jo r t b 2 Jo b 2 n ^ Jl 



r T M dt 



D3n<TM 



D 3n + 2 (if 

D 3 



nb\ 



< irR 2 E (mmb[D 3n , D 3n+2 \y 

D3n<TM 



<nRHt 3 E e-^ 2 ^<nRHt" E exp [ - 3(g) tog(l + n) 
(by Lemma 10, on an event of probability larger than 1 — e ) 



Q£ q 

< 7Ti? 2 6 2 



E(i 

neN 



r 3 ^ < 6^-3 



(for 6 large enough), o 



We establish then the crucial control on a s /b s . 



Lemma 13 Almost surely, if r. is bounded, then a s /b s converges as s — > oo , to some 
random limit f el. Moreover for any e > and any M > 3R/2 there exists a lower 
bound b(a, R, M, e) such that for b > b(a, R, M, e) we have : 



P( sup 



0<S<T M U S 



Oo 

bo 



< 



■ l«o I 
bo 



+ 1) e 2 ) > 1 - e . 



Proof a) Let us prove the second assertion first. 

Set p :=mf{s > 0\b s <bl'h^ a2s }, and r] : = inf {s > 



6o 



> e 



Fix < e < 1/3 and M > 3R/2 . Recall that we set r M : = inf{s | r s = M}, and that, 

for bo > b(M,e) , the event i?i := {/3 > tm} has probability at least 1 — e by Lemma 

. , ,_, ( f T M dt r T M Rdt l2f -2"i , , , -i. 

7, and that the event E 2 :— \ + / 5" < o ^ has probability at least 1 — e 

Wo bi Jo r t bf > 

by Lemma 12. Now, the equation governing a s /b s writes for some Brownian motion W 
and for any s > : 



a -± = a Ji + w 

b s b 



a' 



7s 6t 6t Jo v r t ^ &£/ J 2 h b t b t 



whence 

Kl 



bo 



< 



Oo 

&0~ 



< W* 



a 2 r a -£ 2 dt + * 2 [ 

Jo bf Jo 



* R 



nb 2 



dt 



+ 



Jo b t b t 



(recall that W* denotes the two-sided maximum process of W : W*(s) := max \W\ ([0, s]) 
so that we have on E 1 n E 2 , for b > b(a, R, M, e) and < s < tm A 77 : 





a 




6j 


&o 





t m dt 

! 



TM _R 



r t 6f 



+ 



2 U ° + } Jo b 2 
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< W*[(^ + l) 2 ^- 2 ] + (igl + 1) &r 2 < 2 (M + 1) 6 2 - 1 < , 2 (M + i) , 
on an event E 3 of probability larger than 1 — e , for 6 large enough. 

Hence we find that for b > b(cr, R, M, e) we have rj > tm on E\ fl E 2 fl E 3 , whence : 



sup 

0<S<T M 



a 



<(^ + l)e 2 )>l-3e. 



b) To prove the first assertion of the statement, let us onsider then, for any k G N*, 
the hitting time at ■= inf{s > | b s = b(a, R, k,2~ k )}, which is almost surely finite by 
Lemma 6. Moreover, since r. is bounded, there exists almost surely some random integer 
k such that t k = oo . Applying the strong Markov property to the above, we get 



£ P( sup 

keW 



0-fcATfc<S<Tfc 



> 4 



-k I \ a <T k \ 



+ 1)) < 3 ]T T k < 



oo 



keW 



showing by the Borel-Cantelli lemma that almost surely, there exists some random integer 



k' such that sup 

o"fcAr fc <s<r fc 

surely sup 

s>a k 



< A~ k (^p^- + 1) for k > k' . Hence we get almost 
< A~ k + 1) for k > k := k V k' . This implies at once that 



sup 

*>0"fc„ 



< (2 + 1) , and then that sup 

CTfc s>a k 



< 4-* (2^ + 1) for k > k , 



thereby showing that a s /b s satisfies the Cauchy criterion, and then converges as s — > oo. 

o 

7) End of description of the second case in Theorem 5. 
Lemma 14 Set £> := limsupr s . Then almost surely, if r. is bounded, a./b. converges 

Q V S 



Proof Lemma 13 insures that a s /b s converges to some I e R as s — > oo, and we can 
find a sequence s n of proper times increasing to oo such that g = lim r Sn and T Sn = 

2 

-2 , u-2\ 



for every n, which by the pseudo-norm relation implies : -f 

Of 



(l-^)(r; 2 + 6; 2 ),so 



that using Lemma 6 we find as wanted : £ 



Q ■ 



o 



We prove next that when the relativistic diffusion becomes eventually captured by a 
neighbourhood of the hole, it tends to stay in some asymptotic random plane of 1R 3 . 



Lemma 15 Almost surely, if r. is bounded, the direction b s /b s of the angular momen- 

— * rj 

turn b s converges in S as s — > oo : the trajectories are asymptotically planar. 
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Proof As we already saw in the third part of the proof of Theorem 3, it is easily deduced 
from Proposition 2 that 



«r)=<«.*£)=-(^)£-(w< fa H* s " 



b s > ~v r ,/ V A ' ' S) u s v 2& 2 ' r . 

Then for any fceM and s > A; we have (for some standard Brownian motion 5^) : 

— * 

b 



h 
h 



sup 



j2 Jfc 



oo ^.2 



,2 



t 



+ 



oo CT 2 2 



2 6? 



Now Lemma 6 insures that / -f dt < e~ CT k/2 for fc larger than some finite random 



ik b 2 t 



k , and the Borel-Cantelli lemma insures that B* k (e CT k l 2 ) < e °" fc / 6 for fc larger than 



some finite random k' . Hence for k > k V k' we get sup 

s>k 



2 



showing that b s /b s satisfies almost surely the Cauchy criterion and thus converges in § . 
o 

8) End of the proof of Theorem 5. 

The following lemma proves the statement (ii) of Theorem 5. 

Lemma 16 Fix e e]0, 1[, and suppose r &]R,3R/2[ and T = 0. Fix a, (3 such that 
R < a < r < (3 < 3R/2 . Then for b large enough, with probability larger than 1 — e 
g = limsupr s belongs to ]a, (3\ . 



r x y 1 — j . The function / is continuous and srtictly increasing 

2 



Proof Set /(r) : 
3R 

from [R, — ] onto [0, 3 J^ fl ]- For b large enough, by the pseudo-norm relation we have 

Oo 

bo 



f(r Q ) x \Jl + rl/bl e ]/(«), /(/?)[, and then by Lemma 13 the event 



E := if (a) < inf 

l J V 7 0<S<T, 



3R/2 



a s 










< sup 




< 




0<S<T 3J j /2 







has probability larger than 1 — e. 

Now the pseudo-norm relation implies that f(r s ) l{r s >R} < 7- for any s > . So that 

(by continuity of f' 1 ) r s does never hit /3 on £, and then E C {t 3J?/2 = oo}. Finally 
by Lemma 14 and by continuity of we have f(a) < f(g) < f(/3) and then a < g < f3 
on E . o 



The statement (i) in Theorem 5 is seen as in the proof of Theorem (2,2). Indeed, we 
observed there that if r > 3R/2 , T > 4 + , and A := inf {s > I T s = 2 + }, 
then l/|T sAA | is a supermartingale. Now this implies, if T is large enough: 

(2 + P(A < oo) < liminf E(|T sAA | _1 l{ A<0O }) < liminf E(|T sAA | _1 ) < To" 1 , 
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and then that P( lim r s = +00) > P(A = 00) > 1 - (2 + ) /T > 1 - e . 

Then the irreducibility of the relativistic diffusion follows at once from the proof of 
Theorem (2,2). With Statements (i) and (ii), this proves that both cases in Theorem 5 
indeed occur with strictly positive probability, from any initial condition. Of course we 
already proved the dichotomy in Lemma 5. 

Finally the descrition of the first case in Theorem 5 follows directly from Lemma 2 and 
from Theorem (2,3). o 

3.7 Proof of Corollary 3 

Let us begin by a lemma which specifies at which times the relativistic diffusion ap- 
proaches the top of its excursions outside the hole. 

Lemma 17 Let us consider the following stopping times, for any n G N and any e > : 

D' n := min{s > D 3n+1 \T S = 0} and D e n := min{s > D 3n+1 \T s <eb s }< D' n . 
Then almost surely, when r. is bounded, we have: lim m 1 = Q , and lim r^e = g £ < g , 

where g £ is the unique solution between and ^ of the equation (1 — R/ g e ) g~ 2 = £ 2 — e 2 . 

Proof This is only an additional precision to the second case in Theorem 5 : applying the 
pseudo-norm relation at time D' n , we get 

(1 - R/r D 'J r- D \ ~ (1 - R/r D , n )(r- D l + b~ D \) = a 2 D ,Jb 2 D , n — £ 2 as n -> 00 , 

so that, by Lemma 8, owing to the function / used in the proof of Lemma 16, we must 
have the unique possibility: lim = g. Similarly, at time D e n we get 

(1 - R/r D .) r- D l ~ (1 - R/r D ,)(r D l + b D \) = ^ - ^ — , i 2 - e 2 as n 00 , 
so that we must have the unique possibility : lim r D e = g e . o 

n — >oo n 

We shall need to control the angular contribution around the top of the excursions 
outside the hole. This is the aim of the following lemma. Notations are as above. 

f D 'n 

Lemma 18 We have lim sup / b s ds < „} 2 (^|) 4 g on the intersection of the 
event {g < ^} and of an event of probability larger than 1 — e . 

Proof Set E := {g < ^}. Fix M = M(e) > ^ such that the event E \ {supr. < M} 
has probability smaller than e/8. Set E x = E 1 (e) := {supr. < M}. By Lemma 6 we 
can find some 7 = 7(5) such that, setting Ei = E^e) '■= {b s > 7e ( ^ 1 ~ £ - )o " s for all s > 0}, 
the event Eq \ E2 has probability smaller than e/8. Consider then the stopping time 
v = u(e) at which b. hits the lower bound b(a, R, M, e/4) of Lemma 10, which is finite on 
E 2 . Applying now Lemma 10 and the strong Markov property (applied at time u), we find 
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some integer k = k(e) such that, setting E 3 = E 3 (e) := f) n>k [D 3n+l > (1+ * )ct2 logi 
and E 4 = E 1 fl E 2 fl E 3 , the event E \ E 4 has probability smaller than e/2 . 

Set f3 n := 6d 3 „ + i • We have obviously 



n 



E d C B = B 



(e) := fl {V^ > 7 exp[5^ylog(n/fc)]}. 

n>k 



Consider also for all n > k : B n = B n (e) := Q {\J~f3n > 7 exp[ 2 ^ 1+ £ £ ^ log(ra/ 

m=fc+l 

and let us proceed now somewhat as in Lemma 9, to control the variations of b. and T. 
between D 3n+ i and D' n . 

Set T] n := inf js > D 3n+ i \ \og(b s / (3 n )\ = log2|, and r := inf{s > | r s = M} (which 

is infinite on E 4 ). For D 3n+ i < s < (D 3n+ i + f3~ 1 ^ 2 ) Arj n At and n large enough, we 
have : 



log(& s /&) 



< aWl 



(l + ^)^l+^<± + ^<log2 on 5 n nF n , 



ft! 



where F n := {aW* n [2^l 2 ] < §}. Whence B n n F n C K > (L> 3 n+i + /9" 1/2 ) A r}. 
Recalling Remark 12, we observe then that for large enough n we have 

D 3n+2 A (D 3n+1 + Pn 1/2 ) A r < D 3 n+i + ^f- < £>3n+i + Pn 1/2 on B n H F n , whence for 
any constant a: B n fl F n C {?7 n > (-D3„+2 + cx//3 n ) A r} for large enough n . 

Then for any deterministic q,j E N*, set A ? j := {(j + 1)2 -9 < ^ — C E , and fix 
a deterministic integer m such that each A™j := fl { sup r s < £> + j 2 _9 | satisfies 

s>-D 3m+ i 

P(A ?;j \ A™) < £J . Set r™ := inf{s > D 3m+1 \r s = 3 f-^} (which is infinite on A™J). 



Consider K e n := inf{s > D e n \ T S > j3 n }, and fix some a > (to be specified below) 
The equation governing T. implies that for D e n < s < r\ n A D' n A A £ A (Z)^ + a/f3 n ) A r ( 
and for n large enough we have : 



in 
id 



T s = eb D c n + aW n \ I* (T 2 + l- f)dt] + ^ f T t dt+ f (r t -f)%dt- f 



£-dt 

2rt 



< 2e(3 n + aW:[((3 2 n + l)a/(3 n ] + 3 -£a < 2ef3 n + (3 3 J 4 + 3 -f- a < f3 n 
on F' n := {aW:[2a(3 n } < ^/ 4 }, so that F' n C {A £ „ > ^ A^A (££ + a/(3 n ) A r™ }. 

The equation of T. implies also (in the same way) that, for large enough n and for 

Dl < s < Vn A D' n A A- A {Dl + a//3 n ) A r™ : 



T s < 2e/3 n + <7W'* 2a A 



3^ 



/? 2 

29 v 3i?^ 4 
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while on F' n , for s = D e n + a/(3 n , taking a := ^j- x (^|) 4 , the right hand side above is 
< for n large enough. Hence we find that 

B„nF n nF;c {D' n A r™. < -D^j + a/Aj, for large enough n . 
Therefore so far we get on B n n F n fl , almost surely for large enough ra : 



/ "6. d« < 2a = ^(M)4. 



Moreover, since the standard Brownian motions W n and are independent from 
the cx-field ^ r D 3n+1 , we have : 

nK I ^ 3 „ +1 ) < e-v^/(^ 2 ) and P((i /)c | ^ D3n+i) < e -v^/(^) . 

Therefore 

£ P(F n c n B) < £ n B B < £ ex P ( - I^r exp[^y log("/*)]) < °° , 

n>fc n>fc n>k 

and the Borel-Cantelli lemma implies that B is almost surely included in liminf F n . This 

n 

is obviously the same for F' n , proving that almost surely for any q, j we have on B fl A™- , 
for large enough n : 

f D 'n i , Qe2\ SRx4 



Finally, since the event E \ { limsup / b s ds < „}\ (^) 4 e\ is the increasing union 
(for q e N*) of the events (J A, \ { limsup / " 6 S cfa < ^ (^#) 4 }, we get 



E \ { lim sup / " b s ds < ^ (f 2 ) 4 e}} < F[E \ B] + sup £ F[A qJ \ A™} < e , 

which concludes the proof, o 

We can now establish the second assertion of Corollary 3. (The first assertion is just 
Lemma 17.) 

Set g(r) :=R-r + £ 2 r 3 , so that g(g) = and g'(r) = (VZ \£\r + l)(y/3 \l\r - 1) . 
We have g < ^ <=>■ g < , and g = ^ g = . So that the root r = g of g is 

double if and only if g — ^ . So that the integral / — converges if and only if 

Jo y/ r #(r) 

£> < ^ , thereby justifying our assumption g < ^ . 



47 



Then the the pseudo-norm relation and Lemmas 6 and 14 imply the convergence of 
' ° + (1 — — ) r~ 2 towards I 2 = (1 — — ) g~ 2 . By the definition of b , this can be written : 



U s ds = sgn(T s ) (r s x[R-r s + {f + e s )r s 3 ] ) 1/2 dr s , 

with lim e s = . Lemma 15 means that we can restrict to the limit plane orthogonal to 
(3 := lim b s /b s e S 2 . Indeed, we have 9 S /U S = (b s /b s ) A 9 S , whence 9 S = U s (3 A9 S + e s ) , 
with £ S _L# S and lim e s = ; so that writing # s = X s 3 + V S with K-L3, we have as 
s -> oo : |K| = |/3 A F s | = |/3 A 0„| = l^/t/, - e* s | — ► 1 , whence A s — ► . 

— » — » — > — * — * 

Moreover 3 AV S = (3 A 9 S and then, setting v s := K/|K| , we have on one hand: 



\9 s -v s \ < IAJ + iiyj 



almost surely as s — > oo , 



and on the other hand : 

V s d9 s -d\ s p V s 



dv s = d 



WJ 



IK 



■d\V,\ = A -5- U 9 ds - -t- d\ s + -S- £ s ds - -5- 



= (/3A^)[/ s d S + ( £ -;-[/3Ai; s ])(/3Ai; s )^d S ( since rf{7 s _L{/3, 



(/3A{T S ) (l + det[^- J,v s ])U s ds. 



\v s \ 

Let us now denote by <p s the angular coordinate of v s in the constant plane (3 1 - ; 
then the preceding equation writes equivalently (using the expression of U s seen at the 
beginning of this proof, and choosing the orientation in the plane (3 1 - in order to have a 
positive sign for the remainder of this proof) : 



d<f) s = (i + det[& /3, v s ])U s ds = (l + det[&- /3, v s ]) 



sgn(T s ) <2r s 



yjr a [R-r s + (£ 2 + e 



sjr 3 s ] 



whence (setting 5 S := det , f3 , v s ) for any n6N and e > : 



/ du* 



(1 + 5 S ) dr s 



D, n+1 s JD Sn+1 s JD Sn+1 ^J rs [R-r s + (I 2 + e 8 ) rf] 





I n 




lim 


/ d9 s - 


/ dv s 


n^oo 







by the 



beginning of this proof, we have henceforth to deal with f s s 

Jd ^ Jr a [R-r s + (I 2 + Es) rf] 

(and analogously with D' n instead of D^). Recall that lim r D e = g e < g by Lemma 17. 
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Now, pushing somewhat further the observation we already used for proving Theorem 
3, we can use the strictly increasing radius r s as alternative "time" variable on the proper 
time interval [D 3n+1 , D' n ] . So we set s = s n (r) , e™ := e Sn(r) , <5™ := 5 Sn(r) on [D 3n+1 , D' n ] , 
and we get almost surely : 

Note that we have = |e s ( r )| — * and = 0(\e s ( r )\) — > 0, uniformly on [_D 3n+1 ,Z^J 
as n — > oo , and that we have for large enough n : r^ < (g £ + f?)/2 < Q ■ 

Hence, using the function g introduced at the beginning of this proof, we see that 

< r < r D e n =>- g(r) > g((g £ + £>)/2) > , and then that dominated convergence holds, for 
any e > , showing that we have almost surely 

<-K , ree dr 



lim / d(f) s = [ 

n ^°° JD Zn+1 Jo yV [i? - r + £2 r 3 



. . , f ee dr re dr 

Since it is clear that / = converges to / = as 

J ° ^r[R~r + £ 2 r 3 ] J o yV [R - r + £ 2 r 3 ] 

rK 

e \ , we are left with the remainder / d<p , which we must control uniformly. Note 

that such control is not obvious at all, since we have a root of the denominator precisely at 
time D' n (at which we are at the top of the excursion), according to the fact that around 
this same time the angular move is much more rapid than the radial one. 

Now we have, for large enough n (using that e < -| =>■ (1 — R/g £ ) g~ 2 > — (f )" 2 =>■ 
^>f^r D e>|): 

1 + 8?) dr _ r D ' n dr 



= r« (l + g)dr < 2 p 

JD* Jr n e . r \R — r 4- (ft 4- ?n\ r 3] Jr n e 



r D* n y J r [R-r+ (£ 2 + e")r 3 ] Jr o- n ^ r [R-r + (I 2 + e")r 3 ] 

Jr n e r Z a Jr n e n JD 



>r D e r Jr D e JD^ 

Therefore, applying Lemma 18, we conclude that we have with probability larger than 

- e : 

rDL 



lim sup / d<p s < „ 12 2 (^j=) A £ , and then 



lim sup [ D ' n d<p s - r dr < -± 2 -mfe+ r dr 

n-oc J D3n+1 Jo ^ r [ R _ r + £2 r 3} 3 « 2 <? ^ J Be yV [i? - r + £ 2 r 3] 

which goes to as e \ , while the left hand side does not depend upon e . This proves 
that indeed we have almost surely : 

fK rK re dr 

lim / d6 s = lim / dd> s = I 

n ^°° JDsn+i JD in+1 Jo 



'D 3n+1 n ^°° JD 3n+1 JO J r [Jl- r + £2 r 3] 
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4 Appendix: Study of timelike and null geodesies 



4.1 Timelike geodesies 

The case a = in Section 2.3 corresponds to geodesies ; precisely it is the case of 
timelike geodesies having speed 1. The equations we obtained in Section 2.3 are here 
simplified as follows : 

a , b constant ; r s = (r s — |it) — — 



2 r A s 2 r 2 



Integrating this last equation (after multiplication by f s ) leads simply, up to some constant, 
to the unit pseudo-norm relation *eo g eo = 1 : 



or equivalently to: \s\ = J 



r s 2 = a 2 -(l-g)(l + §), 
r(s) r d r 



r(0) v / a 2 r 2_( 1 _|)( r 2 + 6 2) 

The equations relative to ip s and to ip s yield easily a real constant k such that : 
. 2 _6 2 __ _^ • k 

Vs ~ a a ■ 2 ' Ys — 2 -2 

r]: 77 sin ifs r z s sin ip s 

They are equivalent to the equation j- s {r 2 s 9 s ) = — b 2 r~ 2 9 S . (This is also a conse- 
quence of Proposition 2 in Section 2.4.1 below, taking there a — .) They are also 
equivalent to the constancy of the angular momentum b . This means in particular that 
9 = —U 2 9 — 2(r/r)9 along geodesies, and therefore that every geodesic is included in 
some plane containing the origin {r = 0}. 

Set u := 1/r and P(u) := (1 — Ru)(l + b 2 u 2 ) . The variation of P as r increases 
from R to 00 is as follows : 

- if b < R\^3 : a 2 — P(u) decreases (as r increases from R to 00) from a 2 to a 2 — 1 ; 

- if b > R\f?> : (as r increases from R to 00) a 2 — P(u) decreases firstly from a 2 to 
a 2 — P(ui) , then increases from a 2 — P(ui) to a 2 — P(u 2 ) , then decreases from a 2 — P(u 2 ) 

to a 2 -I, where ^ > m := 1+Vl f R R2 " - > u 2 := l ^ l ~f 2b - > 0. 

Note that for 6 > Ry/S P(ui) = § + ^f(& 2 - 3i? 2 )(l + Vl - 3i? 2 6- 2 ) G ]§, +00 [ , 
and P(u 2 ) = l + ^(b 2 - 3R 2 )(1 - y/1 - 3i? 2 ^ 2 ) e ]|, min{l, P(ui)}[ . 

Besides, , = = is integrable near a simple root of a 2 — P(l/r) , but not 

near a double root. 

Hence, using the relation r 2 = a 2 — P{u) (which does not allow any value for ro), we 
get the following classification of timelike geodesies. (See [L-L], $ 100, problems 1,2, for a 
partial resolution, and [M-T-W], section 25.5, for a more explicit one). 
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Note that we restrict first to the case of S . The extension to the full Schwarzschild 
space S is then easy : see Remark 13 below. 

1. P monotone 

Case 1.1 : b < R^/3 and \a\ > 1 : a 2 — P has no root, so that (s i— > r s ) runs, in one 
direction or in the other one, an increasing trajectory from R to +oo , slowing down (in 
the increasing case) till the limit speed y/a 2 — 1 . In the particular case \a\ = 1 , we have 
r s ~ (9Rs 2 /4y/ 3 as s -> oo . 

Case 1.2 : b < Ry/3 and a 2 < 1 (but b ^ R\/3 or a 2 7^ 8/9) : a 2 - P has a simple 
root, so that (s 1— > r s ) runs a bounded trajectory, which has 2 ends at i?, and increases 
firstly then decreases, with a unique maximum at Rq e]R, 00 [ . r > Rq is impossible. 

Case 1.3 : b = R\/3 and a 2 = P(ui) = 8/9: a 2 — P has a triple root at = 3R, 
so that we have here either a geodesic which runs monotonically (during an infinite time) 
the interval [R, 3R[ , in one direction or in the other one, or a geodesic included in a circle 
(centred at 0) of radius 3R. Indeed, it is easily verified (looking at if) that such circular 
geodesies correspond precisely to multiple roots of a 2 — P . r > 3R is impossible. 

2. P non-monotone 

Case 2.1 : b > R^/3 and a 2 > max{P(wi), 1} , a 2 > P(u\) : a 2 — P has no root, so 
that r does not vanish and we are brought back to the case 1.1, except that there is an 
acceleration phase on the interval r e ^[ . 

Case 2.2.1 : b > R^3 and a 2 < P(u 2 ) : a 2 -Phas a unique (simple) root, r vanishes 
at r = R e]R, ^-[ , and we are brought back to the case 1.2. 

Case 2.2.2 : b > Ry/3 and P(ui) < a 2 < 1 : a 2 — P has a unique (simple) root, r 
vanishes at r = _R 2 e]^, 00 [ , and we are brought back to the case 1.2, alternatively with 
acceleration and slackening phases, the unique maximum being here at R2 ■ 

Case 2.3 : b > R\/3 and a 2 = P(ui) > 1 : a 2 — P has a double root at {r = ^-} (and 
no other root), so that the trajectory (s 1— > r s ) needs an infinite time to reach this level. 
Such geodesies run monotonically, in one direction or in the other one, either the interval 
[R, ^-[, or the interval ]^-,oo] . There are again here also geodesies included in a circle 
(centred at 0) of radius 

Case 2.4 : b > R^3 and 1 < a 2 < P(ui) : a 2 — P has a two (simple) roots, f 
vanishes at r = R E]R, — f and at r = i?i el — , — f ; we are brought back to the case 1.2 
if r e [i?, i?o] ; ^0 €]Ro,Ri[ is impossible; and if r e [i?i,oo], then (s 1— > r s ) runs an 
unbounded trajectory the 2 ends of which are at 00, with a unique minimum at R± . In this 
last subcase we have r s ~ |s| V^ 2 — 1 (and even r s = |s|Va 2 — l + log(l + |s|) +r + o(l) ) 
as s — > 00 if a 2 > 1, and r s ~ (9i?s 2 /4) 1 / 3 if a 2 = 1. Such (projection of) geodesic runs 
approximately a parabola or a branch of hyperbola. 

Case 2.5.1 : b > R\/3 and a 2 = P(ui) < 1 : a 2 — P has a double root at {r = ^-} and 
a simple root at r = R 2 > ^ ■ Hence the trajectory (s 1— > r s ) needs an infinite time to 
reach the level {r = ^-}, so that there are on one hand (projection of) geodesies which run 
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monotonically the interval [R, ^- [ , in one direction or in the other one, and on the other 
hand (projection of) geodesies which run (during an infinite time) from ^ to ^ via a 
unique maximum at R 2 . Again, there are also here geodesies included in a circle (centred 
at 0) of radius r > R 2 is impossible. 

Case 2.5.2 : b > R\/3 and a 2 = P(u 2 ) : a 2 — P has a double root at {r = ^} and a 
simple root at r = R < ^ . Hence we are brought back to the case 1.2 if tq e [-R, i?o], 
and we have also centred circular geodesies of radius r G]-Ro, ^M^, oo[ is impossible. 

Case 2.6 : 6 > R^/3 and P(m 2 ) < a 2 < min{P(M!), 1} : a 2 — P has 3 distinct simple 
roots, so that r vanishes at r = i? G]i?, ^-[ , at r = i?! G]^-, ^] and at r = R 2 G [^, oo[ ; 
we are brought back to the case 1.2 if r G [R, Ro] ; if ?"o G [i?i,i?2], then (s i— > r s ) 
oscillates periodically, increasing from i?i to i?2 then decreasing from R 2 to R±, running 
approximately an ellipse. r <e]R , Ri]U]R 2 , oo[ is impossible. 

Note that the set of limiting radii such that there exists a geodesic which winds 
asymptotically (either from inside or from outside) around a circle (centred at 0) of radius 
— equals the interval l3i?/2, 3R] . Note also that the set of radii — or — of circles (centred 
at 0) which contain geodesies equals the interval ]3i?/2,oo[. 

This gives a geometrical intrinsic meaning to the particular radii 3R and 3R/2 . 



Remark 13 The extension to the full space S of the preceding classification and de- 
scription of timelike geodesies in the restricted space So is more or less straightforward. 
Indeed the results of Theorem 3 are clearly valid for geodesies as well, with the major 
simplification that a s and b s are constant when a — 0. 

Otherwise, there exist timelike geodesies included in the cylinder {r = R} : looking 
at the above and taking a = 0<|A;|<6,we find easily such solutions, which satisfy 

(p s = Arccos^l — k 2 /b 2 sin[±6(s — so)/R 2 ]j . This completes for the space S the picture 
of timelike geodesies we have just drawn above for the strict Schwarzschild space So . 

Note finally the following intrinsic characterisation of the radius R in the space S : R 
is the minimal radius which can be reached by a timelike geodesic which does not hit the 
singularity. 



4.2 Null geodesies 

For these null geodesies, or light rays, the unit pseudo-norm relation is replaced by 

a 2 -f 2 = (l-R/r) r~ 2 . 

Note that the proper time does not make sense any more, so that the new "time" -parameter 
or abscissa A makes sense only up to an afline transform (A i — > gA + q') , and the constant 
parameters a and b do not make sense both anymore, but only their quotient a := a/b. 
This unique "impact parameter" a of the null geodesic is of course a constant of the 
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geodesic. As for the timelike geodesies, the equations relative to 6 = (</?, ip) show that 
they are planar, so that by means of a trivial change of axis we may consider that (p = ir/2 . 

Thus every null geodesic is determined by the equations (in its own plane, the derivatives 
being relative to the abcissa A) (see also [M-T-W], page 674) : 

a constant ; r 2 + (1 — R/r) r~ 2 = a 2 ; i — ct/(l — R/r) ; ip — r~ 2 . 

Eliminating the abcissa A , we get : 

dip = ± r [a 2 - (1 - R/r)r- 2 ]- 1/2 dr ; dt = -ar(l- R/r)- 1 [a 2 - (1 - i?/r)r" 2 ]" 1/2 dr . 

As r increases from to 3R/2, [a 2 — (1 — R/r)r~ 2 ] decreases from infinity to a 2 — 27^2 > 
and then increases from a 2 — ^2 to a 2 as r increases from 3/2/2 to infinity. Therefore, 
owing to their projection on the coordinate r , we find three cases for the null geodesies : 

- Case : \a\ — 3 J^ R '■ The null geodesic can be either included in {r = 3/2/2}, or 
it can be asymptotic to {r = 3/2/2} , either growing strictly from r = to r = 3/2/2, or 
growing strictly from r = 3/2/2 to infinity (or in the reverse direction). 

- Case 1 : \a\ > 3 J 3R '■ The null geodesic runs from to infinity (with an asymptotic 
velocity : r\ ~ aX), or in the reverse direction. 

- Case 2 : |a| < 3 J 3R '■ The null geodesic can either be reminiscent from a parabola, 
with a minimal radius g' > 3/2/2 , or indefinitely oscillate between r = and a maximal 
radius g G ]/2, 3/2/2 [. 

This last sort of null geodesies, which are recurrent at the singularity r = , is the 
most interesting for us here, since the confined diffusion trajectories of Section 2.4.4 seem 
to have their shape asymptotic to the shape of one such geodesic. Indeed, considering 
the impact parameter a such that |a| < 3 J 3R , the maximal radius g < 3/2/2 solves 
a 2 £> 2 + R/g = 1 , so that a stands for t = lim s a s /b s in Theorem (5, 2) of Section 2.4.4, 
and the angular deviation \1/ during each increase from r = to r = g (or decrease from 
r — g to r = ) has exactly the expression found for confined diffusion paths in Corollary 
3 (and Remark 11) of Section 2.4.4. 



geodesic has its graph invariant under the symmetry (r, t) 1— > (r, —t) . Indeed this is clear 
in Sq , and remains true everywhere by analytic continuation. In the Kruskal-Szekeres 
coordinates, this means invariance under the symmetry (u, v) 1— > (u, —v) . The gluing we 
defined at the singularity r = implies at once that null geodesies must also be symmetric 
with respect to {u = v = 0}. 

As a conclusion, we see that the confined null geodesies run indefinitely a closed analytic 
curve, which in the coordinate plane (r, ip) appears as a figure eight (symmetrical and 
centred at the origin), and in the Kruskal-Szekeres coordinate plane (u, v) looks like a pair 
of round brackets (symmetrical with respect to the coordinate axes and joining the two 
branches {r = 0}). 



Moreover, 




shows that each confined null 
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